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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  may  be  converted  to  SI  (metric) 


units  as  follows: 

Multiply 

By 

To  Obtain 

inches 

2.54 

centimetres 

pounds  (force)  per  square  inch 

6894.757 

pascals 

v.vSwv: 
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SPRING  STIFFNESSES  FOR  BEAM-COLUMN  ANALYSIS 
OF  SOIL-STRUCTURE  INTERACTION  PROBLEMS 

CHAPTER  I 

INTRODUCTION 

Structures  of  Interest 

Design  engineers  have  long  been  concerned  with  the  interaction  be¬ 
tween  soils  and  structures.  The  design  of  slabs  on  grade  and  U-frame 
structures  are  controlled  by  the  soil  pressure  acting  on  them.  These 
structures  must  be  designed  to  account  for  the  interaction  between  the 
structure  and  soil.  The  problem  of  determining  the  relationship  between 
the  structural  displacements  and  the  soil  pressures  is  referred  to  as  a 
soi 1 -structure  interaction  problem. 

Potential  Methods 

The  soil -structure  interaction  problem,  in  general,  is  a  nonlinear, 
three-dimensional  problem.  The  load-deformation  curves  for  all  soils  are 
nonlinear  and  are  affected  by  many  factors  such  as  changes  in  pore  water 
pressures  or  cyclic  loading.  Many  structures  have  changes  in  geometry 
which  can  be  analyzed  only  by  considering  a  three-dimensional  behavior. 
Three-dimensional  finite  element  analyses  are  the  only  types  of  analyses 
which  can  address  these  complicated  problems.  When  both  the  three-dimen¬ 
sional  geometry  and  nonlinear  soil  behavior  are  considered,  time  required 
for  preparation  of  data  for  a  finite  element  model  and  computer  time  costs 
become  prohibitive.  For  these  reasons  three-dimensional ,  nonlinear  finite 
element  analyses  are  impractical  for  preliminary  designs. 
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Because  the  cross  sections  of  slabs  on  grade  and  U-frame  structures 
are  relatively  constant  along  the  length,  analysis  for  preliminary  design 
can  be  made  using  a  two-dimensional  plane  strain  slice  of  the  structure. 
Two-dimensional  f ini te  element  analyses  including  nonlinear  soil  behav¬ 
ior  have  been  performed;  however,  these  analyses  are  still  impractical 
for  preliminary  design  because  of  both  cost  and  time. 

For  many  types  of  soi  1 -structure  systems,  the  structure  maybemodel- 
ed  as  an  assemblage  of  beam-coiumn  elements  in  contact  with  the  soi). 

With  this  simplification,  the  major  problem  is  to  represent  the  pressures 
at  the  soi 1 -structure  interface  as  nonlinear  functions  of  the  interface 
displacements.  Winkler's  hypothesis  has  been  used  for  analysis  of  hori¬ 
zontal  beam-column/soil  systems  (2).  However,  Winkler's  hypothesis  does 
not  account  for  the  two-dimensional  behavior  of  the  soil.  The  solution 
of  a  prismatic  beam-column  beam  on  a  Winkler  foundation  loaded  with  a 
uniform  load  yields  a  constant  displacement  across  the  beam.  The  beam 
from  this  model  has  no  rotations  or  moments  and  results  in  a  simple  one¬ 
dimensional  problem.  A  complete  two-dimensional  analysis  of  the  beam-on- 
foundation  indicates  that  displacements  as  well  as  moments  are  not  con¬ 
stant. 

Although  the  so i 1 -structure  interaction  behavior  is  generally  non¬ 
linear,  for  preliminary  design  or  for  short-term  service  loads,  the  soil- 
structure  interaction  behavior  may  be  assumed  to  be  a  linear  relationship 
between  displacements  and  pressures.  Although  much  experimental  data 
have  been  accumulated  from  plate  bearing  tests  to  produce  "coefficients 
of  subgrade  reactions,"  these  coefficients  are  usable  only  with  Winkler's 
hypothesis. 
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Statement  of  the  Problem 


The  objective  of  this  study  was  to  provide  a  simple  means  of  combin 
ing  a  two-dimensional  foundation  with  a  beam-column  structure  which  can 
be  used  for  preliminary  design  of  slabs  on  grade  or  the  base  slab  of  a 
U-frame  structure.  The  study  examines  the  characteristics  of  the  soil- 
structure  interface  and  develops  a  procedure  for  coupling  the  two-dimen¬ 
sional  foundation  to  a  beam-column  model.  The  problem  then  becomes  one 
of  understanding  the  changes  of  the  character i st ics  of  the  model's  stiff 
ness  terms  and  being  able  to  reproduce  these  characteristics  in  a  simple 
and  accurate  procedure. 
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CHAPTER  II 

PREVIOUS  WORK 

The  solution  of  soi 1 -structure  problems  begins  in  1776,  with  the 
work  of  Coulomb  (1).  Coulomb's  theory  provides  a  means  of  calculating 
earth  pressures  against  retaining  structures.  His  theory  assumes  that 
the  structure  is  free  to  move  a  sufficient  amount  to  either  produce  full 
active  or  passive  earth  pressures.  Consequently,  this  theory  is  useful 
for  calculating  only  maximum  pressures.  This  early  work  is  presently  be¬ 
ing  used  for  limit  state  analysis  for  design  of  earth  retaining  struc¬ 
tures.  Thus  it  is  valid  for  predicting  failure  loads,  but  gives  no  in¬ 
formation  about  displacements  or  stresses  between  extremes. 

Winkler  (2),  in  1867,  developed  a  soil  model  for  evaluating  forces 
on  the  structure  which  are  dependent  on  the  displacements  of  the  struc¬ 
ture.  Winkler  assumed  that  the  stress  at  any  point  on  the  surface  of  the 

i 

soil  is  directly  proportional  to  the  structural  displacement  only  at  that 
particular  point.  This  hypothesis  has  been  used  to  aid  in  analyzing  many 
soi 1 -structure  interaction  problems  but  fails  to  account  for  continuity 
of  the  subgrade. 

Subsequently,  a  procedure  was  developed  by  Biot  (3)  for  analysis  of 
an  infinite  beam  on  an  elastic  two-  or  three-dimensional  foundation.  This 
procedure  was  based  on  Winkler's  hypothesis  and  had  no  provisions  for 
plates  or  finite  length  beams  and  considered  only  concentrated  loads. 

Vesic  (5)  extended  Biot's  work  to  include  finite  length  beams  and 
loading  by  a  couple.  Vesic  (4)  also  studied  the  validity  of  Winkler's 
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(i.e.,  the  nodes  on  the  interface)  through  the  middle  nodes  are  plotted. 
These  figures  demonstrate  that  all  relative  terms  of  the  interface  stiff¬ 
ness  matrices  are  very  close  in  magnitude  except  for  the  first  two  (i.e., 
the  end)  rows.  These  figures  also  demonstrate  that  the  terms  change  only 
slightly  for  different  numbers  of  interface  nodes.  These  plots  were  for 
a  foundation  with  a  Poisson's  ratio  of  0.3.  The  same  characteristics  for 
other  values  of  Poisson  were  observed. 

Figures  8,  9,  10,  and  11  are  plots  of  the  major  diagonal  terms  of 
the  four  interface  stiffness  matrices.  These  plots  show  little  change  in 
the  major  diagonal  terms  after  the  first  two  columns  and  little  effect 
from  change  of  number  of  interface  nodes. 

Figure  12  is  a  plot  of  the  centerline  row  of  the  four  i  nterface  sti  ff- 
ness  matrices.  Figure  13  is  a  plot  of  the  first  column  of  the  four  i  nter¬ 
face  stiffness  matrices.  Except  for  the  11-node  matrix,  these  plots  in¬ 
dicate  that  the  magnitude  of  the  stiffness  terms  are  dependent  on  rela¬ 
tive  position  from  the  main  diagonal  and  not  dependent  on  number  of  nodes 
of  the  interface.  These  similarities  of  the  curves  suggest  that  an  inter 
face  stiffness  matrix  for  a  two-dimens ional  foundation  could  be  generated 
without  resorting  to  a  finite  element  solution. 

The  following  is  a  summary  of  the  characteristics  of  the  curves 
shown  in  Figures  4  through  13: 

1.  OFF  diagonal  terms  approach  infinity  as  x  -*■  0  (i.e.,  as  x  ap¬ 
proaches  the  major  diagonal). 

2.  OFF  diagonal  terms  approach  zero  as  x  °°. 

This  suggests  that  the  magnitude  of  the  stiffness  terms  could  be  calcu- 


Poisson's  ratio  for  soils  varies  from  near  zero  to  near  one-half. 

To  examine  the  influence  of  Poisson's  ratio,  interface  matrices  were  de¬ 
veloped  for  Poisson's  ratio  ranging  from  0.1  to  0.4.  Although  magni¬ 
tudes  of  elements  of  stiffness  matrices  depend  on  the  value  of  Poisson's 
ratio,  characteristics  of  matrices  to  be  discussed  later  are  the  same. 
Hence,  only  matrices  for  Poisson's  ratios  of  0.2  and  0.3  are  presented 
and  discussed  in  detail.  To  determine  the  effect  of  the  number  of  nodes 
on  the  interface,  complete  interface  matrices  were  evaluated  for  11,  21, 
31,  and  41  interface  nodes. 

Method  of  Extracting  Interface  Matrix  K| | 

By  definition,  any  element  K.  ^  of  a  stiffness  matrix  is  the  force 
corresponding  to  the  degree  of  freedom  (DOF)  j_  due  to  a  unit  displace¬ 
ment  of  DOF  j_with  all  other  displacements  equal  to  zero.  Because  the 
finite  element  model  of  the  foundation  is  symmetric,  it  is  only  required 
that  the  interface  matrix  be  developed  for  nodes  on  and  to  one  side  of 
the  axis  of  symmetry.  The  interface  matrix  was  developed  by  sequentially 
applying  a  unit  vertical  displacement  at  one  node  with  displacements  of 
other  nodes  on  the  interface  equal  to  zero.  The  reactions  generated  at 
the  restrained  interface  nodes  comprise  one  column  (and  row)  of  the  de- 
si  red  matrix. 

A 

Characteristics  of  K| | 

To  study  K( j ,  interface  matrices  were  generated  for  11,  21,  31,  and 
41  nodes  on  the  interface.  The  off  diagonal  terms  for  each  row  of  each 
matrix  were  plotted  as  shown  in  Figures  4,  5.6,  and  7.  The  terms  below 
the  major  diagonal  term  to  the  end  of  each  column  of  the  first  rows 


Equation  (3.2b)  may  be  written  as: 


or 
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3.3 


3.4a 


3.4b 

3.4c 


where  Kj  j  is  the  interface  stiffness  matrix  which  contains  all  character¬ 
istics  of  the  two-dimensional  foundation. 

The  solution  of  Equation  (3-4)  is  inexpensive  since  the  number  of 

simultaneous  equations  is  the  same  as  for  the  solution  of  the  beam-column 

* 

Once  the  interface  stiffness  matrix,  K(  ( ,  has  been  determined,  the  beam- 

column  stiffness  matrix,  K^,  and  load  array,  Fg,  can  be  changed.  This 

allows  the  solution  of  any  load  case  or  changes  in  the  beam  without  re- 
* 

solving  for  K| j . 


Factors  Affecting  Interface  Matrix  K| | 


Because  the  foundation  is  assumed  to  be  linearly  elastic,  homogene- 

A 

ous,  and  isotropic,  the  interface  matrix  K j  j  is  directly  proportional  to 
the  modulus  of  elasticity  of  the  foundation  material.  Consequently,  in¬ 
terface  matrices  need  be  developed  only  for  a  unit  value  of  foundation 
modulus  of  elasticity.  Other  factors  which  i nfl uence  the  interface  stiff¬ 
ness  matrix  are  the  number  of  nodes  on  the  interface  and  Poisson's  ratio 


for  the  foundation. 


Figure  3.  Window  at  Interface  of  FEM  Foundation  Model 
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incompatible  mode  to  eliminate  parasitic  shear.  The  triangles  are  three 
node  constant  strain  elements.  Figure  3  shows  a  window  display  of  the 
grid  below  the  structure.  For  a  beam  with  a  length  of  20  feet*,  the  grid 
models  an  area  with  a  height  of  3*5  times  and  a  width  of  7  times  the 
length  of  the  beam. 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  (metric) 
units  is  presented  on  page  4. 


Figure  2.  FEM  Foundation  Model  for  20-Foot  Interface 


7 


investigation  that  the  soil  can  be  represented  as  a  homogeneous  isotrop¬ 
ic  semi-infinite  half  space.  Because  preliminary  designs  are  frequently 
based  on  short-term  service  loads  which  result  in  small  displacements, 
it  is  further  assumed  that  the  soil  is  linearly  elastic.  Under  these  as¬ 
sumptions  the  soil  foundation  is  modeled  initially  as  an  assemblage  of 
two-dimensional  finite  elements. 

A  final  simplification  is  made  at  the  interface  between  the  beam- 
column  and  the  soil  foundation.  Enforcement  of  compatibility  of  nodal 
displacements  on  the  interface  would  result  in  shear  stresses  in  the  soil 
along  the  interface  and  a  corresponding  axial  stress  resultant  in  the 
beam-column.  This  effect  would  represent  frictional  interaction  between 
the  base  of  the  structure  and  the  foundation.  The  magnitude  of  this  fric¬ 
tion  resistance  is  dependent  on  roughness  of  the  structure  as  well  as  the 
cohesion  and  internal  friction  of  the  soil.  The  friction  effect  could 
be  accounted  for  by  introduction  of  an  interface  element  (21)  between 
the  structure  and  foundation.  However,  preliminary  studies  indicated 
that  the  interface  friction  had  negligible  influence  on  the  response  of 
the  system.  Consequently,  the  friction  effect  is  neglected  and  compati¬ 
bility  of  vertical  displacements  only  is  enforced  at  the  interface,  which 
allows  the  interface  nodes  to  be  placed  at  the  longitudinal  axis  of  the 
beam-column  rather  than  at  the  base  of  the  beam  with  no  loss  in  accuracy. 

Finite  Element  Model 

Figure  1  displays  the  interface  area  between  the  beam-column  beam 
and  finite  element  foundation;  the  entire  finite  element  foundation  model 
is  shown  in  Figure  2.  The  grid  contains  2139  nodes  and  2124  elements. 

The  quadrilaterals  are  four  node  Isoparametric  elements  using  the 


CHAPTER  I  I  I 


TWO-DIMENSIONAL  FINITE  ELEMENT  ANALYSIS 
OF  THE  STRUCTURE  SOIL  SYSTEM 

Objective 

The  objective  of  this  study  is  to  provide  a  simple  means  of  combin¬ 
ing  a  two-dimensional  foundation  with  a  beam-column  structure  which  can 
be  used  for  preliminary  design  purposes.  To  accomplish  this  goal,  a  pro¬ 
cedure  to  produce  an  interface  model  which  approximates  the  results  ob¬ 
tained  from  a  finite  element  solution  was  developed. 

Beam-Column/Finite  Element  Method 


Assumptions 

As  stated  previously,  attention  is  limited  to  those  structures  such 
as  slabs-on-grade  or  U-frame  base  slabs  for  which  a  two-dimensional  slice 
of  the  soi 1 -structure  system  adequately  represents  the  behavior  for  pre¬ 
liminary  design.  It  is  further  assumed  that  the  dimensions  of  the  struc¬ 
ture  allow  a  model  composed  of  beam-column  elements  to  be  used  to  repre¬ 
sent  this  part  of  the  system. 

Inherent  in  the  two-dimensional  model  of  the  soi 1 -structure  system 
is  the  assumption  that  the  soil  foundation  is  in  a  state  of  plane  strain. 
Although  the  extent  of  the  soil  foundation  away  from  the  structure  will 
influence  the  soi 1 -structure  interaction,  it  is  assumed  for  this 


dimensional  finite  element  model  of  the  foundation  indicates  that  the  mo¬ 
ments,  rotations,  and  relative  displacements  are  accurately  predicted  by 
the  continuum  model.  However,  the  absolute  displacements  are  substantial 
ly  in  error. 

The  finite  element  method  has  proven  to  be  a  valid  tool  for  analyz¬ 
ing  soi 1 -structure  interaction  problems.  Clough  and  Duncan  (19)  used  the 
finite  element  method  to  analyze  sheet  pile  walls  and  U-frame  structures 
that  were  a  part  of  the  Port  Allen  and  Old  River  Locks.  The  procedure  ac 
counts  for  nonlinear  stress-strain  behavior  of  the  soil  medium  and  uses 
an  incremental  loading  technique  to  simulate  a  construction  sequence.  The 
favorable  comparison  of  analytical  and  experimental  results  for  the  Port 
Allen  and  Old  River  Locks  demonstrated  the  usefulness  of  this  method  in 
the  soil-structure  interaction  area  (19,  20).  However,  the  use  of  this 
method  requires  the  tedious  generation  of  complicated  grids  for  modeling 
the  problem,  large  amounts  of  computer  time,  and  many  hours  of  interpret¬ 
ing  the  results  of  analysis.  Therefore,  the  high  cost  of  this  analysis 
technique  makes  it  impractical  for  use  by  most  designers. 


the  Fi lonenko-Borodich  model.  The  Pasternak  (12,  14)  model  uses  a  shear 
layer  to  Interact  between  the  Winkler  springs  and  the  structure.  Vlasov 
and  Leont'ev  (15,  16)  developed  a  two-parameter  foundation  model  from  the 
continuum  point  of  view. 

All  two-parameter  models  can,  in  general,  be  represented  by  the  fol¬ 
lowing  equation*: 

q(x)  -  K  w (x)  -0  — 4X) 
s  dxz 

where 

q(x)  =  surface  pressure; 

w(x)  =  surface  displacements; 

K  =  modulus  of  subgrade  reaction;  and 
=  second  soil  parameter. 

The  two  parameters,  K  and  D,  are  difficult  to  determine.  Presently,  these 
parameters  are  being  determined  from  extensive  field  tests.  Although  the 
parameter  K  has  the  same  name  as  the  K  for  Winkler's  model,  i.e.,  modulus 
of  subgrade  reaction,  the  values  for  K  for  a  Winkler  model  cannot  be  used 
for  a  two-parameter  model.  This  foundation  model,  when  combined  with  a 
structural  model,  leads  to  a  very  complex  mathematical  problem.  Because 
of  the  difficulties  in  calculations  and  in  determining  the  two  parameters, 
this  is  presently  not  a  practical  tool  for  preliminary  design. 

Elastic  continuum  models  allow  for  a  continuous  description  of  the 
soil  beneath  the  structure  as  well  as  displacements  of  the  soil  away  from 
the  structure.  Flamant  (17)  and  Cheung  (18)  calculated  a  flexibility  ma¬ 
trix  for  the  interface.  The  inverse  of  this  flexibility  matrix  yields  a 
stiffness  matrix  for  the  foundation  which  could  be  combined  with  a  beam- 
column  stiffness  matrix.  A  comDarison  of  this  approach  with  a  two- 


Symbols  are  listed  and  defined  in  the  Notation  (Appendix  N) 


hypothesis  and  concluded  that  the  hypothesis  is  valid  for  infinite  beams 
resting  on  a  semi-infinite  elastic  subgrade.  Vesic  stated  that  the  use 
of  Winkler's  hypothesis  leads  to  an  overestimation  of  bending  moments  and 
an  underestimation  of  contact  pressures  and  displacements.  Vesic  also 
noted  that  conventional  analysis  using  Winkler's  hypothesis  could  be  used 
for  finite  beams  and  recommended  procedures  for  beams  with  different  length 
characteristics. 

Reese  and  Matlock  (6)  used  the  Winkler  hypothesis  with  a  nonlinear 
foundation  for  analysis  of  laterally  loaded  piles.  While  the  foundation 
was  a  nonlinear  Winkler  model,  the  structure  was  modeled  with  beam-column 
elements.  The  finite  difference  technique  was  used  to  formulate  the  beam- 
column  equations. 

Haliburton  (7)  extended  the  procedure  used  by  Matlock  et  al.  (8)  for 
analyzing  flexible  retaining  structures.  Haliburton  discussed  nonlinear 
soil  response,  anchor  and  brace  supports,  and  procedures  for  developing 
nonlinear  soil  load-deflection  curves. 

Oawkins  (9)  combined  the  nonlinear  Winkler  foundation  with  a  stiff¬ 
ness  method  approach  to  the  beam-column  model.  This  work  also  allowed 
the  use  of  linearly  varying  distributed  springs  as  well  as  point  linear/ 
nonl i near  springs. 

Two-parameter  models  were  developed  in  order  to  overcome  deficien¬ 
cies  inherent  in  the  Winkler  hypothesis  in  modeling  the  continuous  behav¬ 
ior  of  the  subgrade.  The  two-parameter  models  improved  on  Winkler's 
model  by  adding  a  continuous  layer  between  the  structure  and  the  Winkler 
springs.  The  Fi lonenko-Borodich  model  (10,  11,  12)  adds  a  membrane  under 
tension  between  the  structure  and  the  Winkler  spring.  The  Hetenyi  (12, 

13)  model  adds  an  elastic  beam  for  two-dimensional  problems  in  place  of 


where  y  is  the  magnitude  of  the  interface  stiffness  term,  and  x  is  the 
number  of  terms  from  the  major  diagonal; 


f(x)  -*•  0  as  x  -*■  0 
f (x)  -  BxC 

where  A,  B,  C  are  functions  of  Poisson's  ratio. 


Equilibrium  Requirement 

it 

Each  row  of  K j  j  is  a  system  of  forces  produced  by  a  unit  displace¬ 
ment  of  a  node  on  the  interface  with  all  other  nodes  of  the  interface 
node  having  no  displacements.  The  forces  produced  on  the  interface  nodes, 
together  with  the  reactions  produced  at  the  boundary  nodes  of  the  com¬ 
plete  finite  element  model,  constitute  a  system  of  forces  in  equilibrium. 
If  the  rows  of  the  interface  matrix  are  approximated,  the  terms  must  con¬ 
tinue  to  satisfy  overall  equilibrium.  Therefore, 


n 


I 

i*>l 


j  *  1  ,  2 ,  3 ,  •  «  • ,  N 


where 


n 

y 


K 


'  J 


number  of  node  on  interface; 
approximated  row;  and 

finite  element  interface  matrix  stiffness  coefficient. 


CHAPTER  IV 


GENERATION  OF  APPROXIMATE  INTERFACE 
AND  COMPARISON  OF  RESULTS 


Introduction 

Because  the  stiffness  matrix  of  the  beam  column  is  independent  of 
the  interface  stiffness  matrix,  any  beam  can  be  used  for  comparison  of 
results.  For  the  following  comparison,  a  beam  12  inches  high  and  1  inch 
thick  will  be  used.  Poisson's  ratio  of  0.3  and  Young's  modulus  of  3x10 
will  be  used  for  the  material  properties  of  the  beam.  The  four  beams 
will  have  length-to-depth  ratios  of  5.  10,  15»  and  20.  The  node  spacing 
for  all  beams  will  be  six  inches  in  order  to  accurately  represent  the 
beam.  Because  each  beam  node  is  connected  to  the  interface,  the  inter¬ 
face  node  spacing  is  also  six  inches.  This  node  spacing  will  allow  for 
accurate  modeling  of  the  soi 1 -structure  interaction  behavior.  The  fol¬ 
lowing  is  a  list  of  the  five  different  load  cases  used  for  comparison: 
Load  Case  1 — Uniform  load  over  entire  beam 
Load  Case  2--Point  load  at  centerline 
Load  Case  3 — Point  load  at  left  quarter  point 
Load  Case  k — Uniform  load  over  center  one-half  of  beam 
Load  Case  5"~Uniform  load  over  left  half  of  beam. 


To  provide  a  basis  for  comparison  with  approximate  interface  matrices, 


displacements  and  rotations  using  the  exact  interface  matrices  were  ob¬ 
tained  for  Poisson's  ratios  of  0.1,  0.2,  0.3»  and  0.4  for  I)  and  21 
node  beams  using  all  load  cases.  The  exact  displacements  for  the  31  and 
41  were  calculated  for  Poisson's  ratios  of  0.2  and  0.3  using  all  five 
load  cases. 

For  reinforcement,  complete  solutions  of  the  beam-soil  system  using 
two-dimensional  finite  elements  for  both  beam  and  foundation  were  obtain¬ 
ed  for  11-  and  21-node  beams  with  a  foundation  having  a  Poisson's  ratio 
of  0.3  and  a  Young's  modulus  of  7000  pounds  per  square  inch  (psi).  The 
7000  psi  represents  a  dense  sand  (22).  The  results  were  compared  with 
solutions  for  a  beam-column  model  of  the  beam  with  a  two-dimensional 
fi ni te  element  foundation.  As  expected,  these  solutions  produced  identi¬ 
cal  results  for  both  models.  Additional  reinforcement  resulted  from  the 
solutions  for  11- ,  21- ,  and  31-node  beams  coupled  to  a  41-node  exact  inter¬ 
face  matrix.  As  expected,  these  solutions  produced  the  same  results  as 
when  each  beam  was  coupled  with  an  interface  matrix  containing  the  same 
number  of  nodes.  An  example  of  a  21  beam  on  a  21-node  foundation  compar¬ 
ed  with  a  21  beam  on  a  4l-node  foundation  for  a  uniform  load  over  the  en¬ 
tire  beam  is  shown  in  Table  1.  A  Poisson's  ratio  of  0.3  was  used  for  the 
foundations  for  this  comparison. 


Expedient  Approach 


Because  the  beam-column  matrices  and  the  foundation  matrices  are  in¬ 
dependent,  it  is  possible  to  produce  a  beam-column  computer  program  which 
stores  complete  interface  matrices  for  various  Poisson's  ratios.  Only 
one  interface  matrix  for  each  Poisson's  ratio  is  required,  provided  the 
matrix  is  sufficiently  large  to  accommodate  a  wide  variety  of  beam 
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TABLE  1 


21-NODE  BEAM  ON  41-NODE  INTERFACE 


Node  Number 
From  Left 

End  of  Beam 

Beam  on  4l- 

Node  Foundation 

Beam  on  21- 

Node  Foundation 

Deflection 

(in.) 

Rotation 

Deflection 

(in.) 

Rotation 

1 

-0.00245166 

-0. 312634x1 0”^ 

-0.00245166 

-0.312637xl0”5 

2 

-0.00247034 

-0.308966x10”^ 

-0.00247034 

-0.308967xl0”5 

3 

-0.00248859 

-0.2981 1 IxlO-5 

-0.00248859 

-0.2981 12xl0”5 

4 

-0.00250596 

-0.279648x10-5 

-0.00250597 

-0.279649x10"5 

5 

-0.00252201 

-0.254003x10 

-0.00252201 

-0.254004x10 

6 

-0.00253631 

-0.221918x10-5 

-0.00253632 

-0.221919xl0"5 

7 

-0.00254853 

-0. 184316x10”'* 

-0.00254853 

-0.184317x10"5 

8 

-0.00255834 

-0.l42222xl0"5 

-0.00255835 

-0. I42223xl0”5 

9 

-0.00256552 

-0.96721 5x1 0“6 

-0.00256553 

-0.967217x10”6 

10 

-0.00256990 

-0. 489344x1 O”6 

-0.00256991 

-0.489345x10"6 

t 

-0.00257137 

0.0 

-0.00257138 

0.0 

lengths.  A  61-node  interface  would  cover  a  majority  of  beam  lengths. 

The  user  of  such  a  program  would  be  required  to  provide  only  the  modulus 
of  elasticity  and  Poisson's  ratio  for  the  beam  and  foundation,  length  of 
the  beam,  moment  of  inertia  of  the  beam,  and  loading  on  the  beam.  This 
seems  feasible  since  only  961  terms  are  needed  to  produce  a  61-node  in¬ 
terface  matrix,  as  the  matrices  are  symmetric  about  both  diagonals.  The 
studies  of  the  interface  matrices  indicated  that  the  magnitude  of  terms 
in  any  row  exhibited  identical  variations  with  position  from  the  major 
diagonal  as  shown  in  Figures  4  through  7.  Consequently,  this  indicates 
the  possibility  of  reproducing  the  interface  behavior  with  fewer  than 
961  terms. 

Generation  of  Matrix  From  Centerline  Row 

Introduction 

This  centerline  row  method  develops  the  entire  interface  matrix  us¬ 
ing  only  the  centerline  row  of  a  61-node  exact  interface  matrix  and  30 
terms  to  force  equilibrium.  This  procedure  generates  an  approximate  in¬ 
terface  matrix  producing  results  sufficiently  accurate  for  preliminary 
designs  for  beams  with  up  to  41  nodes.  Beams  with  more  than  41  nodes 
will  have  larger  errors  because  the  beam  action  will  be  influenced  by  the 
outer  rows.  As  seen  in  Appendix  A,  the  changes  in  any  interface  matrix 
from  one  row  to  the  next  are  small,  except  for  the  outside  few  rows. 

This  end  effect  of  the  outside  row  is  difficult  to  predict  and  is  not 
modeled  by  this  procedure.  This  method  models  the  matrix  by  using  the 
centerline  row  and  forcing  equilibrium  of  all  rows. 


As  previously  stated,  each  row  of  the  interface  matrix  represents  a 
set  of  reactions  resulting  from  a  single  interface  node  being  displaced 
a  unit  value.  The  sum  of  all  vertical  reactions  is  zero,  i.e.,  the  sum 
of  vertical  reactions  of  the  interface  nodes  and  of  the  boundary  of  the 
finite  element  grid.  The  sum  of  all  horizontal  forces  is  also  zero. 

Thus,  the  problem  is  in  a  State  of  equilibrium. 

The  approximate  matrix  must  also  satisfy  equilibrium.  If  the  approx¬ 
imated  interface  matrix  does  not  satisfy  equilibrium,  the  matrix  will  not 
produce  reasonable  results.  For  example,  a  10  percent  change  of  the 
smallest  off  diagonal  term  of  a  21-node  interface  matrix  connected  to  a 
21-node  beam  will  result  in  a  3  percent  error  of  all  displacements.  This 
indicates  that  the  results  produced  from  these  interface  matrices  are 
very  sensitive  to  the  requirement  for  satisfying  of  equilibrium.  A  ma¬ 
trix  which  does  not  satisfy  equilibrium  results  in  a  change  in  the  inter¬ 
nal  strain  energy  of  the  foundation  and  misrepresents  the  true  system  of 
forces.  Thus,  each  approximate  interface  matrix  must  be  generated  such 
that  each  row  satisfies  equilibrium. 

Each  row  can  be  forced  to  satisfy  equilibrium  by  making  the  sum  of 
each  approximate  row  equal  the  sum  of  the  corresponding  row  of  the  exact 
interface  matrix.  The  tables  in  Appendix  B  give  the  sums  of  terms  on  each 
row  of  the  interface  stiffness  matrices  listed  in  Appendix  A.  The  sum  of 
each  row  was  divided  by  the  sum  of  terms  of  the  centerline  row  to  produce 
factors  which  can  be  used  to  force  equilibrium  of  generated  matrices  by 
knowing  the  sum  of  the  centerline  row  for  a  given  Poisson's  ratio.  Figure 
I1*  displays  the  ratios  of  the  sums  of  the  rows  plotted  for  each  row  for  a 
Poisson's  ratio  of  0.2.  Figure  15  displays  the  ratios  of  the  first  term 


plotted  against  the  number  of  nodes  on  the  interface.  This  illustrates 
that  the  relationship  (i.e.,  ratio)  of  the  sum  of  the  first  row  with  the 
sum  of  the  centerline  row  changes  very  little  with  increased  nodes  on 
the  interface.  These  figures  suggest  that  equilibrium  factors  could  be 
extrapolated  from  these  data.  However,  extrapolation  of  these  data  would 
produce  results  with  errors  greater  than  those  which  will  be  reported 
later  in  this  study. 

1>1  Node  Foundation 

Table  2  lists  the  results  of  using  the  previously  described  beam 
with  21  nodes  on  a  4l-node  interface  matrix,  [E],  for  Poisson's  ratio  of 
0.3  which  was  generated  by  the  centerline  row  method.  The  maximum  dis¬ 
placement  due  to  a  uniform  load  is  1.1  percent  in  error,  while  the  maxi¬ 
mum  moment  is  2.0  percent  in  error.  These  results  indicate  that  the  cen¬ 
terline  row  method  of  producing  an  approximate  interface  method  can  be 
used  for  preliminary  design.  The  following  is  the  centerline  row  method 
for  generating  the  interface  matrix: 

1.  Place  the  centerline  row,  which  depends  on  the  desired  Poisson's 
ratio,  in  every  row  starting  at  the  major  diagonal  (Figure  16).  The  place¬ 
ment  of  the  centerline  row  results  in  a  banded  matrix  with  a  half  band¬ 
width  equal  to  half  the  number  of  terms  of  any  row  plus  one,  since  all 
interface  matrices  have  an  odd  number  of  nodes.  The  odd  number  of  nodes 
for  each  interface  forces  a  node  to  be  at  the  axis  of  symmetry. 

2.  Force  symmetry  about  the  major  diagonal. 

3.  Satisfy  equilibrium  on  the  first  row  above  the  centerline  row 
(row  20).  This  row  contains  one  unknown  term,  E£q  ,  where 
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E20,<,1  '  ExaCt  E20,  i  ‘  ,1,  A«>rC,x-  E20 ,  i 

The  exact  matrix  [E ]  has  the  same  Poisson's  ratio  as  the  approximate  ma¬ 
trix. 

4.  Force  symmetry  about  the  minor  diagonal: 

El,22  =  E20,41 

5.  Satisfy  equilibrium  on  row  19.  The  row  contains  two  unknowns, 

E 1 9 , 40  and  E 1 9 , 4 1 : 

41  39 

Exact  E 1 9 , i  "  ApprOX-  E 1 9 ,  i  =  E 1 9 , 40  +  E 1 9 , 4 1 

Because  values  at  the  end  of  each  row  in  the  matrix  become  essentially 
constant  (Figures  4  through  7),  it  is  assumed  that  Ejg  =  Ejg  ^ . 

6.  Force  symmetry  about  minor  diagonal: 

E1,23=EI9,41  and  E2,23  =  E19,40 

7.  The  process  is  continued  for  each  successive  row  above  the  cen¬ 
terline  row.  Each  row  contains  one  or  more  unknown  terms.  When  more 
than  one  unknown  is  involved,  all  the  involved  terms  are  assigned  the 
same  value  for  that  row. 

8.  Each  term  which  has  been  previously  determined  by  forcing  sym¬ 
metry  about  the  minor  diagonal  must  not  be  changed.  For  example,  row  two 
has  19  unknown  terms,  after  placement  of  centerline  row  and  forcing  sym¬ 
metry  about  the  major  diagonals: 

E£  2  through  E2  ^ 2  =  Centerline  row 
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design  and  analysis,  it  is  preferable  over  the  second  procedure,  which 
approximates  the  off  diagonal  terms  with  an  exponential  equation.  The 
second  procedure  produces  acceptable  results  but  has  at  this  point  few 
additional  benefits  other  than  requiring  slightly  fewer  numerical  values. 

Recommendat ions 

Both  procedures  developed  an  approximate  semi- i nf i ni te  half  space. 
The  effect  of  finite  depths  on  this  soi 1 -structure  interaction  problem 
was  not  investigated.  The  effect  for  shallow  depth  is  a  common  problem 
for  foundations  and  needs  to  be  investigated.  The  procedure  which  uses 
an  equation  to  produce  the  off  diagonal  terms  would  be  more  beneficial 
if  the  terms  of  the  equation  could  be  related  directly  to  depth,  Pois¬ 
son's  ratio,  and  the  modulus  of  elasticity.  This  goal  could  possibly  be 
achieved  by  investigating  other  equations  with  more  coefficients.  Fur¬ 
ther  studies  based  on  the  theory  of  elasticity  solutions,  such  as  Fla- 
mant's,  could  lead  to  finding  the  correct  relationship  between  the  coef¬ 
ficients  of  the  equation  and  the  physical  properties  of  the  foundation. 


CHAPTER  V 


SUMMARY,  CONCLUSION,  AND  RECOMMENDATIONS 
Summary 

The  intent  of  this  research  was  to  investigate  the  soi 1 -structure 
interaction  for  a  beam-column  structure  on  a  continuous  linearly  elastic 
subgrade,  and  to  retain  the  simplicity  of  the  beam-column  analysis  while 
providing  a  foundation  model  which  includes  the  two-dimensional  behavior 
of  the  foundation.  The  goal  was  to  develop  a  simple  method  to  produce  a 
two-dimensional  foundation  behavior  which  can  be  used  for  preliminary  de 
sign  and  analysis  of  slabs  on  grade  and  U-frame  structures. 

This  research  showed  that  the  foundation  could  be  investigated  and 
formed  independently  of  the  structural  beam-column  matrix.  Also,  it  was 
shown  that  the  interface  matrix,  which  represents  a  semi-infinite  half 
space,  could  be  generated  from  a  few  terms.  The  terms  can  be  determined 
from  known  values  of  Poisson's  ratio  and  the  modulus  of  elasticity  of 
the  foundation.  Two  procedures  to  generate  a  soi 1 -structure  interaction 
stiffness  matrix  were  illustrated. 

Conclusion 

The  first  procedure,  which  generates  the  interface  stiffness  matrix 
from  the  centerline  row,  produces  the  better  results  and  has  One  less  ap¬ 
proximation.  Since  this  procedure  satisfies  the  goal  of  producing  an 
easily  constructed,  two-dimensional  foundation  useful  for  preliminary 
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This  procedure  illustrates  that  an  interface  stiffness  matrix  can 
be  developed  from  36  terms  ([four  coefficients.  A,  B,  C,  and  D] ,  31  sums 
of  rows,  and  the  major  diagonal  term  of  centerline  row).  These  36  terms 
are  dependent  only  on  Poisson's  ratio  and  the  modulus  of  elasticity. 

This  procedure  produces  sufficiently  accurate  results  for  preliminary  de 
signs  with  25  fewer  terms  by  assuming  an  equation  for  the  off  diagonal 
terms.  This  procedure  could  easily  be  included  into  a  beam-column  analy 
sis  computer  program  for  analysis  and  design  of  slabs  on  grade  and  U- 
frame  structures.  The  procedure  has  the  same  limitation  as  the  previous 
procedure  of  not  accounting  for  foundations  with  specific  depths  instead 
of  the  results  of  a  half  space  solution. 


SUMMARY  OF  RESULTS 


row  procedure,  the  selection  of  the  major  diagonal  term  depends  only  on 
the  desired  Poisson's  ratio. 

_gvC 

2.  Use  the  equation  (y  -  Ae  +  D)  to  calculate  the  first  30  off 
diagonal  terms.  The  values  for  the  coefficients  A,  B,  C,  and  D  for  Pois¬ 
son's  ratios  of  0.1,  0.2,  0.3»  and  0.4  are  given  in  Table  4. 

3.  Force  symmetry  about  the  major  diagonal. 

4.  Force  every  row  of  the  generated  stiffness  matrix  to  satisfy 
equilibrium  according  to  the  sums  in  Appendix  B.  The  sums  chosen  from 
this  appendix  depend  on  the  value  of  Poisson's  ratio.  The  procedure  was 
described  in  detail  for  the  41 -node  matrix. 

5.  Construct  the  matrix  to  be  symmetrical  about  the  minor  diagonal. 
The  procedure  for  forcing  symmetry  was  also  previously  described  for  the 
41-node  interface. 

6.  Multiply  the  entire  matrix  by  the  desired  modulus  of  elasticity. 

Appendices  K  thfough  M  give  the  results  of  the  five  load  cases  for 

a  Poisson's  ratio  of  0.3  and  a  Young's  modulus  of  7000  psi  for  21- ,  31- , 
and  41-node  beams  using  the  procedure  described  above.  Table  5  gives  a 
summary  of  the  errors  in  maximum  deflection,  moment,  and  shear  using 
a  Poisson's  ratio  of  0.3  and  the  five  load  cases. 

These  results  demonstrate  the  same  characteristics  as  the  results 
from  the  centerline  row  method.  The  errors  increase  as  the  number  of 
beam  nodes  (length  of  beam)  increase.  This  increase  in  error  for  this 
procedure  is  a  result  of  the  end  effects  from  interface  matrices  as  was 
the  reason  for  the  centerline  row  method.  However,  the  errors  in  the 
results  for  the  coefficient  method  are  greater  than  those  from  the  previ¬ 
ous  results,  since  the  procedure  has  an  added  assumption  of  the  approxi¬ 
mation  of  the  off  diagonal  terms. 
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makes  one  additional  approximation  in  that  procedure.  An  equation  is 
used  to  calculate  the  off  diagonal  terms  of  the  centerline  row  of  the 
interface  stiffness  matrix. 

Figures  4  through  7  display  the  characteristic  shape  of  the  rows  of 
each  interface  stiffness  matrix.  After  examination  of  several  equations, 
the  one  found  to  best  fit  the  data  is  of  the  form: 

-BXC 

^(stiffness)  +  ® 

where  X  is  the  number  of  nodes  from  major  diagonal,  and  B,  C,  and  D  are 
constant.  An  example  of  this  fit  is  shown  in  Figure  17. 

The  above  equation  cannot  be  fitted  to  the  stiffness  values  by  the 
conventional  least  squares  method.  Therefore,  a  sequential  search  method, 
described  in  Appendix  J,  was  used  to  find  the  coefficients  A,  B,  C,  and 
0  of  the  equation.  The  fit  is  accomplished  by  letting  the  sequential 
method,  a  simplex  method  by  Nelder  and  Mead  (23),  vary  A,  B,  C,  and  D 
until  the  sum  of  the  squares  of  differences  between  the  FEM  stiffness 
values  and  the  calculated  stiffness  values  becomes  a  minimum. 


Procedure 

This  procedure,  the  coefficient  method,  is  identical  to  the  center- 

line  row  method  except  for  the  approximation  of  the  off  diagonal  terms 

of  the  centerline  row.  By  knowing  the  general  shape  of  the  variation  of 

_gj(C 

the  off  diagonal  terms  (y  =  Ae  +  D) ,  the  following  procedure  can  be 
used  to  develop  an  approximate  6l-node  interface  stiffness  matrix: 

1.  Make  all  major  diagonal  terms  a  single  value,  i.e.,  the  first 
term  of  the  centerline  row  listed  in  Appendix  C.  As  with  the  centerline 
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constant  displacement,  rigid  body  rotation  (Appendices  D,  E,  G,  and  H.) 

The  second  reason  is  that  the  depth  of  foundation  to  the  length  of  the 
beam  becomes  smaller  with  larger  beams.  The  depth-to-length  ratio  for  the 
41-node  beam  is  3.5,  which  means  the  interface  is  still  an  approximation 
for  an  elastic  half  space.  When  the  depth-to-length  ratio  becomes  small¬ 
er,  the  assumption  of  an  elastic  half  space  is  not  valid  and  the  inter¬ 
face  will  begin  to  be  influenced  by  the  boundary  of  the  finite  element 
grid. 

Conclusion 

This  procedure  illustrates  that  an  approximate  interface  stiffness 
matrix  can  be  developed  from  the  centerline  row  of  the  exact  matrix  and 
the  ratios  required  to  satisfy  equilibrium.  Therefore,  a  6l-node  inter¬ 
face  matrix  can  be  produced  from  61  terms  (31  term  centerline  row  of  ma¬ 
trix  and  30  ratios).  The  procedure  is  more  acceptable  for  inclusion  wi th 
a  beam-column  program  than  the  previous  procedure  because  of  the  fewer 
terms  required  to  produce  the  interface  matrices.  The  procedure  still 
has  the  limitation  of  not  accounting  for  foundations  with  specific  depths 
rather  than  the  half  space  as  used  herein. 

Generation  of  Matrix  With  Approxi¬ 
mated  Centerline  Row 

Introduct ion 

It  has  been  shown  that  the  6l-node  interface  stiffness  matrix  gener¬ 
ated  from  the  centerline  row  can  be  used  for  beams  up  to  41  nodes  and 
thus  eliminates  the  need  for  developing  interface  stiffness  matrices  for 
each  beam  with  a  different  number  of  nodes.  The  following  procedure 


foundation.  Also,  the  percentage  in  the  maximum  moment  has  decreased 
from  2.0  percent  error  for  a  41 -node  centerline  row  interface  to  0.84 
percent  error  for  a  61 -node  centerline  row  interface.  This  reduction  in 
error  is  due  to  the  effect  of  the  outside  row  of  the  approximate  inter¬ 
face  being  further  away  from  the  end  of  the  beam. 

The  shear  diagrams  in  Appendices  D  through  I  have  a  stair-stepped 
shape  for  sections  of  the  beam  where  loads  have  not  been  applied,  and 
have  a  saw-toothed  shape  for  sections  which  have  been  loaded  with  a  uni¬ 
formly  distributed  load.  These  discontinuities  occur  at  each  beam  node 
and  the  magnitude  of  these  discontinuities  are  the  reactions  produced  by 
the  foundation  at  that  particular  node.  The  elastic  half  space  founda¬ 
tion  has  been  discretized  and  is  represented  by  the  interface  matrix 
which  acts  as  a  series  of  interconnected  linear  springs  attached  to  each 
beam  node.  These  discontinuities  will  become  smaller  as  the  node  spac¬ 
ing  of  the  foundation  becomes  smaller. 

As  seen  in  Table  3,  this  method  produces  good  results  for  beams 
with  up  to  41  nodes.  Also,  this  table  illustrates  that  as  the  number  of 
nodes  increases,  so  does  the  percentage  error.  The  6l-node  centerline 
row  interface  is  limited  to  a  maximum  of  approximately  41  nodes  for  two 
reasons.  The  first  reason  is  demonstrated  in  Table  3,  which  shows 
greater  errors  with  the  4l  nodes,  particularly  with  unsymmetric  loads 
near  the  edge  of  the  beam  (Load  Case  5).  Appendices  F  and  I  illustrate 
that  for  41-node  beams  the  deformed  shape  for  Load  Case  5  (uniform  load 
over  left  half  of  the  beam)  is  different  from  the  correct  deformed  shape 
The  difference  in  this  case  is  not  constant;  some  disp’ snts  are 
above  the  correct  value  while  others  are  below.  The  difference  in  the 
deformed  shape  for  smaller  beams  for  other  load  cases  are  in  error  by  a 
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Therefore,  E ^  23  through  ^  are  unknown  after  forcing  symmetry  about 
the  major  diagonal.  Terms  E£  23  through  E^  ^  w‘ 11  be  determined  from 
forcing  symmetry  about  minor  diagonal.  E^  23  =  Eig  i*0’  E2  2k  ~  E|8  40’ 

E2,25  ’  E17,40  •  •  •  E2,39  “  E3,W  Therefore,  only  terms  E^  end  Ej,,, 
are  unknown  when  approximate  row  2  is  forced  to  satisfy  equilibrium. 

9.  Multiply  the  entire  matrix  by  the  desired  modulus  of  elasticity 
This  41-node  approximated  matrix  can  be  used  only  with  beams  of  a  limit¬ 
ed  number  of  nodes.  Beams  with  more  than  21  nodes  will  have  greater 
errors  than  the  2  percent  error  in  moment  reported  above.  This  is  a  re¬ 
sult  of  improper  modeling  of  the  outside  rows  of  the  interface  matrix. 

61- Node  Foundation 

By  using  a  6l-node  foundation,  the  troublesome  interface  end  effects 
are  moved  further  away  from  the  end  of  the  beam.  Appendix  C  gives  the 
centerline  row  for  Poisson's  ratio  of  0.1,  0.2,  0.3,  and  0.4.  Appen¬ 
dices  D  through  J  give  the  results  of  the  five  load  cases  for  Young's 

modulus  of  7000  psi  with  Poisson's  ratios  of  0.2  and  0.3  for  21-,  31-,  ai 

> 

41 -node  beams.  Table  3  gives  a  summary  of  the  errors  of  maximum  de¬ 
flection,  moment,  and  shear  for  two  different  Poisson's  ratios  and  the 
five  load  cases.  The  6l-node  interface  matrix  generated  for  these  re¬ 
sults  was  produced  using  the  steps  outlined  for  the  41-node  generated 
i nterface. 

Appendix  D  presents  the  results  of  the  21-node  beam  on  a  6l-node 
centerline  row  interface  matrix  for  Poisson's  ratio  of  0.3.  As  seen  in 
Table  3,  the  percentage  error  of  the  maximum  displacement  for  a  uniform 
load  over  the  entire  beam  is  0-91  percent.  This  percentage  error  is  bet 
ter  than  the  1.1  percent  error  reported  for  the  same  beam  on  a  4l-node 


APPENDIX  A 


INTERFACE  STIFFNESS  MATRICES 


Poisson  Ratio  0.2 

-0.4687744  0.2014B56  0.0801580  0.0289753  0.0200524  0.0138527  0.0107073  0.0088199  0.0077943  0.0082193 
0.0113046 

-0.5551036  0.1472307  0.0678919  0.0205690  0.0143566  0.0095844  0.0073999  0.0061676  0.0061966  0.0082193 
-0.5686461  0.1624127  0.0644176  0.0183880  0.0127684  0.0084279  0.0065749  0.0061674  0.0077963 
-0.5703274  0.1613199  0.0639394  0.0179566  0.0125569  0.0084279  0.0073999  0.0088199 
-0.5710193  0.1409287  0.0437420  0.0179566  0.0127684  0.0095844  ..0107073 
-0.5711841  0.1609287  0.0639394  0.0183880  0.0143566  0.0138527 


Poi sson  Ratio  0. 3 

-0.4964217  0.2130380  0.0852976  0.0305544  0.0207889  0.0144729  0.0111426  0.0091479  0.0080547  0.0084601 
0.0115481 

-0.5877406  0.1766444  0.0723771  0.0218514  0.0148399  0.0100300  0.0077026  0.0064041  0.0064124  0.0084601 
-0.6022470  0.1715333  0.0689623  0.019SSS4  0.0131726  0.0088129  0.0068396  0.0064041  0.0030547 
-0.4039966  0.1704024  0.0682524  0.0191044  0.0129469  0.0088129  0.0077026  0.0091479 
-0.6046983  0.1700019  0.0680524  0.0191044  0.0131724  0.0100300  0.0111426 
-0.4048700  0.1700019  0.0682524  0.0195554  0.0148399  0.0144729 


A1-- 11-Node  Interface 


-0.4497093  0.2008391  0.0795837  0.0283694  0.0193791  0.0130671  0.0097499  0.0075756  0.0041079  0.0050546 

0.0042791  0.0036886  0.0032324  0.0028769  0.002601]  0.0023934  0.0022496  0.0021827  0.0022039  0.0026336 

0.0039857 

-0.5555509  0.1668329  0.0674720  0.0201019  0.0138109  0.0089184  0.0065326  0,0049883  0.0039730  0.0032580 

0.0027381  0.0023493  0.0020540  0.0018289  0.0014404  0.0015423  0.0014807  0.0014810  0.0017531  0.0026336 

-0.5490201  0.1620386  0.0642007  0.0179004  0.0121723  0.0076494  0.0055139  0.0041496  0.0032686  0.0026570 

0.0022187  0.0018960  0.0016557  0.0014780  0.0013526  0.0012811  0.0012663  0.0014810  0.0022039 

-0.5707233  0.1608783  0.0434220  0.0173227  0.0117269  0.0072940  0.0052236  0.0039084  0.0030659  0.0024854 

0.0020734  0.0017744  0.0015549  0.0014034  0.0013114  0.0012811  0.0014807  0.0021827 

-0.5715127  0.1603494  0.0430309  0.0170223  0.0114881  0.0071003  0.0050641  0.0037761  0.0029560  0.0023951 

0.0020014  0.0017210  0.0015247  0.0014034  0.0013526  0.0015423  0.0022496 

-0.5718656  0.1600897  0.0628326  0.0168660  0.0113630  0.0069989  0.0049820  0.0037104  0.0029051  0.0023590 

0.0019816  0.0017210  0.0015569  0.0014780  0.0016604  0.0023934 

-0.5720557  0.1599457  0.0627204  0.0167776  0.0112930  0.0069441  0.0049406  0.0036819  0.0028897  0.0023590 

0.0020014  0.0017744  0.0016557  0.0018289  0.0026011 

-0.5721636  0.1598631  0.0626567  0.0167290  0.0112571  0.0069200  0.0049280  0.0036819  0.0029051  0.0023951 

0.0020736  0.0018960  0.0020540  0.0028769 

-0.5722254  0.1598171  0.0626236  0.0167071  0.0112461  0.0069200  0.0049406  0.0037104  0.0029560  0.0024854 

0.0022187  0.0023493  0.0032324 

-0.5722576  0.1597964  0.0626134  0.0167071  0.0112571  0.0069441  0.0049820  0.0037761  0.0030659  0.0026570 
0.0027381  0.0036886 

-0.5722676  0.1597964  0.0626236  0.0167290  0.0112930  0.0069989  0.0050641  0.0039084  0.0032686  0.0032580 
0.0042791 


A2--21-Node  Interface  With  Poisson  Ratio  of  0.2 


C~7 


I 


-0.497520? 

0.0044394 

0.0039241 

-0.5881771 

0.0028450 

-0.4025943 

0.002304? 

-0.4043910 

0.0021549 

-0.4051954 

0.0020821 

-0.4055421 

0.0020417 

-0.4037390 

0.0020821 

-0.4056704 

0.002154? 

-0.4059344 

0.002304? 

-0.4059477 

0.0028450 

-0.4059781 

0.0044394 


0.2124117  0.0847381 
0.0038143  0.0033337 


0.0299431  0.0201257  0.0134959  0.0101914 
0.0029547  0.0024424  0.0024387  0.0022801 


0.1742541 

0.0024343 

0.1711423 

0.0019440 

0.149939? 

0.0018387 

0.1494159 

0.0017850 

0.1491444 

0.0017850 

0.1489974 

0.0018387 

0.1489120 

0.0019440 

0.1488444 

0.0024343 

0.1488431 

0.0038163 

0.1688431 


0.0719429 

0.0021214 

0.0685466 

0.0017094 

0.0477317 

0.0014077 

0.0473303 

0.0015760 

0.0471246 

0.0016077 

0.0670067 

0.0017096 

0.0669429 

0.0021214 

0.0669087 

0.0033337 

0.0668981 


0.0213879  0. 
0.0018821  0. 
0.0190671  0, 
0.0013203  0 
0.0184640  0, 
0.0014437  0 
0.0181557  0. 
0.0014437  0 
0.0179940  0 
0.0015203  0 
0.0179026  0 
0.0018821  0 
0.0178324  0 
0.0029347 
0.0178300  0 


0142961  0. 
0017017  0. 
0123731  0, 
0013853  0 
0121084  0. 
0013421  0 
0118639  0. 
0013833  0, 
0117343  0 
0017017  0, 
0116621  0 
,0026626 
0116251  0, 


0093637 

0015731 

0080301 

0013047 

0076599 

0013047 

0074417 

0013731 

0073570 

0024387 

0073006 


0.0068347 

0.0015007 

0.0057659 

0.0012829 

0.0054639 

0.0013007 

0.0053011 

0.0022801 

0.0052164 


0.0079113 

0.0021971 

0.0052166 

0.0014923 

0.0043383 

0.0014923 

0.004087? 

0.0021971 

0.0039330 


0.0063639 

0.0022043 

0.0041467 

0.0017337 

0.0034106 

0.0022043 

0.0032003 


0.0032390 

0.0026160 

0.0033936 

0.0026160 

0.002744? 


0072737  0.0051609 


.0072757  0.005173? 


0.0178300  0.0116251  0.0073006  0.0052164 
0.0178524  0.0116621  0.0073370  0.0053011 


A3 — 21-Node  Interface  With  Poisson  Ratio  of  0.3 


-0.4698703 

0.0041450 

0.0013053 

0.001883? 

-0.5556184 

0.0026390 

0.0008053 

-0.5690657 

0.0021250 

0.0006271 

-0.3707664 

0.0019673 

0.0005760 

-0.5715571 

0.0018736 

0.0005497 

-0.5719137 

0.0018193 

0.0005424 

-0.3721103 

0.0017847 

0.0005497 

-0.3722271 

0.0017627 

0.0005760 

-0.5723013 

0.001748? 

0.0004271 

-0.5723306 

0.0017413 

0.0008053 

-0.3723840 

0.0017390 

0.0013033 

-0.5724070 

0.0017413 

-0.3724227 

0.001748? 

-0.3724327 

0.0017427 

-0.3724384 

0.0017847 

-0.3724403 

0.0018193 


0.2007330  0.0794986  0.0282869  0.0192956  0.0129801  0.0096576 
0.0035370  0.0030584  0.0026749  0.0023419  0.0021033  0.0018871 
0.0012090  0.0011276  0.0010397  0.0010053  0.0009654  0.0009429 


0.0074760  0.0059991 
0.0017049  0.0013501 
0.0009471  0.0009863 


0.1667774 

0.0022357 

0.0007470 

0.1619941 

0.0017871 

0.0005831 

0.1608344 

0.0016484 

0.0003383 

0.1603033 

0.001565? 

0.0005181 

0.1600384 

0.001317? 

0.0003181 

0.139886? 

0.0014876 

0.0005383 

0.1597936 

0.0014687 

0.0005831 

0.1597331 

0.0014374 

0.0007470 

0.1396924 

0.0014523 

0.0012090 

0.1596446 

0.0014323 

0.1396434 

0.0014374 

0.1594327 

0.00146B7 

0.139424? 

0.0014874 

0.1596211 

0.001517? 

0.1596211 

0.001365? 


0.0674180 

0.0019219 

0.0006987 

0.064133? 

0.0015266 

0.0003479 

0.0633763 

0.0014039 

0.0005100 

0.0629817 

0.0013306 

0.0004974 

0.0627771 

0.0012881 

0.0003100 

0.0626560 

0.001261? 

0.000547? 

0.062379? 

0.0012437 

0.0006987 

0.0625296 

0.0012370 

0.0011276 

0.0624956 

0.0012343 


0.0200473 

0.0016724 

0.0006600 

0.0178536 

0.0013214 

0.0005216 

0.0172741 

0.0012121 

0.0004921 

0.0149691 

0.001146? 

0.0004921 

0.0168056 

0.0011093 

0.0005216 

0.0167064 

0.0010844 

0.0006600 

0.0166430 

0.0010733 

0.0010597 

0.0166007 

0.0010671 

0.0165720 

0.0010671 


.0137540 

.0014707 

.0006314 

.0121226 

.0011569 

.0005054 

.0116746 

.0010590 

.0004889 

.0114303 

.0010007 

.0005034 

.0112961 

.0009677 

.0006314 

.0112134 

>.0009481 

.0010053 

.011159? 

.0009377 


0.0088581 

0.0013053 

0.0006146 

0.0075957 

0.0010229 

0.0005041 

0.0072367 

0.0009351 

0.0005041 

0.0070363 

0.0008830 

0.0006146 

0.0069241 

0.0008540 

0.0009654 

0.0068541 

0.0008377 

0.0068084 

0.0008303 


0.0064676 

0.0011681 

0.0006154 

0.0054551 

0.0009126 

0.0005217 

0.0051603 

0.0008336 

0.0006154 

0.0049926 

0.0007871 

0.0009429 

0.0048974 

0.0007620 

0.0048374 

0.0007490 

0.0047981 

0.0007450 


0.0049173 

0.0010531 

0.0006391 

0.0040849 

0.0008211 

0.0006391 

0.0038379 

0.0007499 

0.0009471 

0.0036951 

0.0007086 


0.0038944 

0.0009564 

0.0007857 

0.0031963 

0.0007446 

0.0009863 

0.0029860 

0.0006804 

0.0028631 

0.0006441 


0.0624723  0.0165524 
0.0012370  0.0010733 


0.0111240  0.006777?  0.0047719 
0.0009344  0.0008303  0.0007490 

0.0110996  0.0067570  0.0047543 
0.0009377  0.0008377  0.0007620 

0.0110831  0.0067433  0.0047430 
0.0009481  0.0008540  0.0007871 


0.0036134  0.0027924 
0.0006873  0.000626? 

0.0035617  0.0027473 
0.0006779  0.0006216 

0.0035277  0.002717? 
0.0006779  0.000626? 

0.0035051  0.0026984 
0.0006873  0.0006441 

0.0034903  0.0026863 
0.0007086  0.0006804 

0.0034813  0.0026794 
0.000749?  0.0007446 


0.0624564 

0.0012457 

0.0624461 

0.0012619 

0.0624404 

0.0012881 

0.0624384 

0.0013304 

0.0624404 

0.001403? 


0.0165393 

0.0010864 

0.0165313 

0.0011093 

0.0165276 

0.001146? 

0.0165276 

0.0012121 

0.0145313 

0.0013214 


0.0110724 

0.0009677 

0.0110666 

0.0010007 

0.0110647 

0.0010590 

0.0110666 

0.001156? 

0.0110726 

0.0014707 


0.0067350 

0.0008830 

0.0067310 

0.0009351 

0.0047310 

0.001022? 

0.0067350 

0.0013053 

0.0067433 

0.0021033 


0.0047367 

0.0008336 

0.0047347 

0.0009126 

0.0047367 

0.0011681 

0.0047430 

0.0018871 

0.0047543 


0.0034771 

0.0008211 

0.0034771 

0.0010531 

0.0034813 

0.001704? 

0.0034903 


0.0026773 

0.0009544 

0.0026794 

0.0015501 

0.0026863 


0.0049364 

0.0014181 

0.0012154 

0.0031704 

0.0008744 

0.0012154 

0.0025753 

0.0006806 

0.002394C 

0.0006230 

0.0022873 

0.0005916 

0.0022254 

0.000578? 

0.0021860 

0.000578? 

0.0021604 

0.0005916 

0.0021440 

0.0006230 

0.0021341 

0.0006806 

0.0021294 

0.0008744 

0.0021294 

0.0014181 

0.0021341 


Ai» — 31-Node  Interface  with  Poisson  Ratio  of  0.2 


-0.4974413  0.3123200  0.0844424  0.0298893  0.0200304 
0.0043144  0.0034744  0.0031707  0.0027444  0.0024344 
0.0013184  0.0012131  0.0011270  0.0010324  0.0009914 
0.0017324 

-0.3882374  0.1742043  0.0719141  0.0213381  0.0142440 
0.0027513  0.0023244  0.0019957  0.0017327  0.0013200 
0.0008147  0.0007520  0.0004994  0.0004344  0.0004234 
-0.4024374  0.1711217  0.0485051  0.0190237  0.0123247 
0.0022130  0.0018393  0.0015850  0.0013489  0.0011934 
0.0004347  0.0005871  0.0005484  0.0005189  0.0004991 
-0.4044309  0.1499191  0.0474890  0.0184184  0.0120589 
0.0020517  0.0017141  0.0014587  0.0012544  0.0010953 
0.0005831  0.0005423  0.0005109  0.0004894  0.0004823 
-0.4052371  0.1493721  0.0472834  0.0181031  0.0118083 
0.0019547  0.0014324  0.0013844  0.0011909  0.0010347 
0.0005574  0.0005224  0.0004989  0.0004894  0.0004991 
-0.4054080  0.1490973  0.0470713  0.0179354  0.0114494 
0.0019010  0.0015834  0.0013413  0.0011527  0.0010031 
0.0005301  0.0005224  0.0003109  0.0005189  0.0004234 
-0.4058114  0.1489404  0.0449440  0.0178331  0.0113840 
0.0018437  0.0015524  0.0013144  0.0011294  0.0009833 
0.0005574  0.0003423  0.0005484  0.0004544  0.0009914 
-0.4039324  0.1488440  0.0448474  0.0177477  0.0115290 
0.0018433  0.0015333  0.0012981  0.0011140  0.0009727 
0. 0005831  0.0005871  0.0004994  0.0010524 
-0.4040091  0.1487814  0.0448134  0.0177241  0.0114920 
0.0018291  0.0015219  0.0012893  0.0011099  0.0009494 
0.0004347  0.0007520  0.0011270 

-0.4040599  0.1487397  0.0447804  0.0174947  0.0114470 
0.0018214  0.0013144  0.0012844  0.0011099  0.0009727 
0.0008147  0.0012131 

-0.4040943  0.1487110  0.0447544  0.0174744  0.0114501 
0.0018190  0.0015144  0.0012893  0.0011140  0.0009833 
0.0013184 

-0.4041180  0.1484914  0.0447403  0.0174411  0.0114393 
0.0018214  0.0013219  0.0012981  0.0011294  0.0010031 
-0.4041340  0.1484783  0.0447299  0.0174329  0.0114331 
0.0018291  0.0013333  0.0013144  0.0011527  0.0010347 
-0.4041444  0.1484703  0.0447240  0.0174490  0.0114311 
0.0018433  0.0015324  0.0013413  0.0011909  0.0010933 
-0.4041303  0.1484444  0.0447220  0.0174490  0.0114331 
0.0018437  0.0013834  0.0013844  0.0012544  0.0011954 
-0.4041521  0.1484444  0.0447240  0.0174529  0.0114393 
0.0019010  0.0014324  0.0014387  0.0013489  0.0015200 


0.0134170  0.0101073  0.0078200  0.0042454  0.0031477 
0.0021439  0.0019334  0.0017429  0.0015790  0.001438' 
0.0009447  0.0009141  0.000*077  0.0009343  0.0011354 

0.0093080  0.0047743  0.0051503  0.0040731  0.0033110 
0.0013453  0.0012003  0.0010737  0.0009740  0.0008889 
0.0004020  0.0005944  0.0004131  0.0007444  0.0011354 
0.0079799  0.0037104  0.0042770  0.0033417  0.0024887 
0.0010343  0.0009379  0.0008410  0.0007400  0.0004919 
0.0004931  0.0005057  0.0004131  0.0009343 
0.0074054  0.0054034  0.0040200  0.0031234  0.0025010 
0.0009444  0.000B573  0.0007487  0.0004950  0.0004337 
0.0004931  0.0005944  0.0009077 

0.0074000  0.0052319  0.0038743  0.0029983  0.0023923 
0.0009123  0.0008109  0.0007274  0.0004590  0.0004027 
0.0004020  0.0009141 

0.0072841  0.0051337  0.0037901  0.0029254  0.0023287 
0.0008830  0.0007854  0.0007041  0.0004414  0.0005900 
0.0009447 

0.0072120  0.0050720  0.0037370  0.0028791  0.0022884 
0.0008444  0.0007724  0.0004944  0.0004343  0.0005900 

0.0071450  0.0050314  0.0037021  0.0028490  0.0022423 
0.0008590  0.0007483  0.0004944  0.0004414  0.0004027 

0.0071334  0.0050044  0.0034790  0.0028293  0.0022454 
0.0008590  0.0007724  0.0007041  0.0004590  0.0004337 

0.0071123  0.0049844  0.0034439  0.0028147  0.0022354 
0.0008444  0.0007854  0.0007274  0.0004950  0.0004919 

0.0070981  0.0049751  0.0034547  0.0028097  0.0022307 
0.0008830  0.0008109  0.0007487  0.0007400  0.0008889 

0.0070894  0.0049487  0.0034503  0.0028074  0.0022307 
0.0009123  0.0008573  0.0008410  0.0009740  0.0014387 
0.0070854  0.0049444  0.0034503  0.0028097  0.0022354 
0.0009444  0.0009379  0.0010787  0.0015790 
0.0070854  0.0049487  0.0034547  0.4028147  0.0022454 
0.0010543  0.0012003  0.0017429 

0.0070894  0.0049751  0.003443V  0.0028293  0.0022423 
0.0013453  0.0019354 

0.0070981  0.0049844  0.0034790  0.0028490  0-0022884 
0.0021439 


A5--31-Node  Interface  With  Poisson  Ratio  of  0.3 


-0.4*00200  o. 20070 jj 

0.0041210  0.0033117 
0.0012'30  0.001137? 
0.0003013  0.00033(7 
0.0000480 

-0.3330373  0.1007014 
0.0020224  0.0022170 
0.0007301  0.0000030 
0.0003000  0.0003420 
-0.3000707  0.1010830 
0.0021104  0.0017711 
0.0003840  0.0003327 
0.0002774  0. 0002040 
-0.3707703  0.1008230 
0.0010310  0.0010317 
0.0003270  0.0004800 
0.0002333  0.0002420 
-0.3713004  0.1002030 
0.0018370  0.0013470 
0.0004027  0.00044(7 
0.0002407  0.0002320 
-0.3710234  0.1000284 
0.0018007  0.0014074 
0.0004717  0.0004200 
0.0002373  0.0002320 
-0.3721200  0.130(700 
0-0017030  0.0014043 
0.0004380  0.0004171 
0.0002407  0.0002420 
-0.3722384  0.1307817 
0.00173(1  0.0014414 
0.0004400  0.0004003 
0.0002333  0.0002040 
-0.3723140  0.1307107 
0.0017201  0.0014231 
0.0004433  0.0004040 
0.0002774  0.0003420 
-0.3723047  0.1300773 
0.0017070  0.0014134 
0.0004401  0.0004023 
0.0003000  0.00033(7 


0.0704730  0.0282007 
0.0030314  0.0020431 
0.0010431  0.0000033 
0.0003303  0.0003001 

0.0074033  0.0200331 
0.0010020  0.0010313 
0.0000300  0.0003800 
0.0003200  0.0003131 
0.0041430  0.0178430 
0.0013003  0.0013024 
0.0004870  0.00044(0 
0.0002333  0.0002437 
0.0033000  0.0172044 
0.0013(30  0.0011020 
0.0004401  0.0004040 
0.0002330  0.0002313 
0.0020710  0.0100300 
0.0013100  0.0011247 
0. 0004100  0.0003774 
0.00022(0  0.0002313 
0.0027007  0.0107040 
0.0012030  0.0010(41 
0.0003031  0.0003010 
0.0002330  0.0002437 
0.0020440  0.0100043 
0.0012337  0.0010373 
0.0003(20  0.0003317 
0.0002333  0.0003131 
0.0023071  0.0100200 
0.0012131  0.00103(7 
0.0003731  0.0003400 
0.0003200  0.0003001 
0.0023134  0.0103(34 
0.0012000  0.0010230 
0.0003714  0.0003433 
0.0003303 

0.0024703  0.0103344 
0.0011000  0.0010101 
0.0003703  0.0003433 


0.0102707  0.0120010 
0.0023204  0. 0020070 
0.0008000  0.0008257 
0.0004(07  0.0004720 

0.0137410  0.008(437 
0.0014470  0.0012707 
O.OOOS413  0.0003017 
0.0003041  0.0002007 
0.0121123  0.0073831 
0.0011330  0.0000007 
0.0004140  0*0003630 
0.0002420  0.0002450 
0.0110443  0.0072200 
0.0010300  0.0000103 
0.0003733  0.0003400 
0.0002330  0.0002450 
0.0114100  0.0070240 
0.0000701  0.0008354 
0.00034(4  0.0003231 
0.0002420  0.0002007 
0.0112844  0,0000117 
0.0000300  0.0008221 
0.0003341  0.0003103 
0.0003041  0.0004720 
o.oirooo  0. 0008403 
0.0000131  0.0008000 
0.0003230  0.0003030 
0.0004(07 

0.0111454  0.0007020 
0.000(0(3  0.0007840 
0.0003200  0.0002007 

0.0111070  0.0007300 
0.0008(04  0.0007741 
0.0003104  0.0402007 

0.0110800  0.0007303 
0.000(781  0.0007007 
0.0003200  0.0003030 


0.00003(0  0.0074300 
0.0018477  0.0010011 
0.0007070  0. 0007100 
0.0004003  0.0004721 

0.0004347  0.0040037 
0.0011307  0.0010214 
0.0004004  0.0004331 
0.0003034  0.0003180 
0.0054440  0.0040730 

o.oooeaoo  o. 0007021 
0.0003304  0.0003327 
0.0002300  0.0003180 
0.0031480  0.003(230 
0.0008050  0.0007184 
0.0003210  0.0003007 
0.0003034  0.0004721 
0.0040804  0.0030821 
0. 0007330  0.0000730 
0.0003011  0.0002820 
0.0004003 

0.0048840  0.0035001 
0.0007237  0.0000434 
0.0002807  0.0002721 

0.0048220  0.0033430 
0.0007050  0.0000271 
0.0002830  0.0002470 

0.0047(11  0.0033003 
0.0000010  0.0000147 
0.0002(20  0.0002070 

0.0047323  0.0034(37 
0. 0004823  0.0000001 
0.0002(30  0.0002721 

0.0047310  0.0034034 
0.0000757  0.0000004 
0.0002807  0.0002(20 


0.00307(0  0.0040147 
0.0013013  0.0013031 
0.0000000  0.0000281 
0.0004044  0.0000114 

0.0030801  0.0031359 
0.0000207  0.0008341 
0.0004071  0.0003820 
0.0003030  0.0000114 
0.0031830  0.0023017 
0.0007110  0.0000431 
0.0003117  0.0002033 
0.0004044 

0.0020720  0.0023700 
0.0000440  0.0003817 
0.0002823  0.0002000 

0.0028401  0.0022720 
9.0000030  0.0003430 
0.0002050  0.0002517 

0.0027700  0.0022084 
0.0003770  0.0005200 
0.0002374  0.0002437 

0.0027207  0*0021007 
0.0003013  0.0003030 
0.0002350  0.0002437 

0.0020077  0.0021383 
0.0095300  0.0004034 
0.0002574  0.0002517 

0.0020730  0.0021181 
0.0005420  0.0004(10 
0.0002030  0.0002000 

0.0020387  0.0021034 
0.0005370  0.0004(40 
0.0002(23  0.0002033 


A6— 41-Node  Interface  With  Poisson  Ratio  of  0.2  (Cont.) 


-0.5724001  0.1500471  0.0024534 
0.0010074  0.0014040  0.0011024 
0.0004301  0.0004023  0.0003714 
0.0005013 

-0.5724257  0.1500253  0.0024340 
0.0010004  0.0013080  0.0011771 
0.0004401  0.0004040  0.0003751 
-0.5724444  0.1500001  0.0 024204 
0.0010850  0.0013040  0.0011730 
0.0004433  0.0004003  0.0003820 
-0.5724583  0.1505070  0.0024000 
0.0010823  0.0013020  0.0011710 
0.0004400  0.0004171  0.0003031 
-0.5724080  0.1505880  0.0024020 
0.0010804  0.0013007  0.0011710 

0.0004580  0.0004200  0.0004100 
-0.5724700  0.1505811  0.0023001 
0.0010707  0.0013007  0.0011710 
0.0004717  0.0004487  0.0004401 
-0.5724824  0.1505703  0.0023020 
0.0010804  0.0013020  0.00U736 
0.0004027  0.0004800  0.0004870 
-0.5724800  0.1503727  0.0023801 
0.0010823  0.0013040  0.0011771 
0.0005270  0.0005327  0.0000300 
-0.5724804  0.1505700  0.0023874 
0.0010850  0.0013080  0.0011824 
0.0003840  0.0000030  0.0010431 
-0.3724011  0.1505004  0.0023870 
0.0010004  0.0014040  0.0011000 
0.0007301  0.0011370 
-0.3724010  0.1505004  0.0023874 
0.0010074  0.0014134  0.0012000 
0.0012430 


0.0105320  0.0110010  0.0007101 
0.0010004  0.0008723  0.0007017 
0.0003400  0.0003250  0.0003103 

0.0105150  0.0110400  0.0007003 
0.0010049  0.0008084  0.0007580 
0.0003317  0.0003341  0.0003231 
0.0105033  0.0110357  0.0000909 
0.0010020  0.0008001  0.0007570 
0.0003010  0.0003484  0.0003400 
0.0104940  0.0110270  0.0000897 
0.0010000  0.0008054  0.0007370 
0.0003774  0.0003733  0.0003830 
0.0104871  0.0110217  0.0000840 
0.0010000  0.0008001  0.0007580 

0.0004040  0.0004140  0.0005017 
0.0104821  0.0110173  0.0000810 
0.0010020  0.0008084  0.0007017 
0.0004480  0.0003413  0.0008237 
0.0104780  0.0110144  0.0000787 
0.0010049  0.0008723  0.0007007 
0.0005800  0.0008900 
0.0104703  0.0110129  0. 0000770 
0.0010094  0.0008781  0.0007741 
0.0009033 

0.0104753  0.0110123  0.0000770 
0.0010101  0.0008804  0.0007849 

0.0104753  0.0110129  0.0000787 
0.0010230  0.00089B3  0.0008000 

0.0104703  0.0110144  0.0000810 
0.0010387  0.0009151  0.000B221 


0.0047103  0.0034519  0.0020407 
0.0000714  0.0005970  0.0003349 
0.0003011  0.0003007  0.0003117 

0.0047050  0.0034419  0.0026377 
0.0006690  0.0005953  0.0005340 
0.0003219  0.0003327  0.0004071 
0.0046967  0.0034344  0.0026313 
0.0006683  0.0005953  0.0005349 
0.0003566  0.0004351  0.0006696 
0.0046904  0.0034290  0.0026266 
0.0006690  0.0005970  0.0005376 
0.0004664  0.0007160 
0.0046860  0.0034253  0.0026236 
0.0006714  0.0006004  0.0003426 

0.0007679 

0.0046830  0.0034230  0.0026217 
0.0006757  0.0006061  0.0005500 

0.0046813  0.0034219  0.0026211 
0.0006823  0.0006147  0.0005613 

0.0046807  0.0034219  0.0026217 
0.0006919  0.0006271  0.0005779 

0.0046813  0.0034230  0.0026236 
0.0007056  0.0006434  0.0006030 

0.0046830  0.0034253  0.0026266 
0.0007257  0.0006730  0.0006446 

0.0046B60  0.0034290  0.0026313 
0.0007539  0.0007184  0.0007119 


0.0020927 

0.0004829 

0.0003826 

0.0020849 

0.0004829 

0.0006281 

0.0020791 

0.0004849 

0.0020753 

0.0004889 


0.0020727 

0.0004954 


0.0020716 

0.0005056 

0.0020716 

0.0003206 

0.0020727 

0.0003436 

0.0020753 

0.0005817 

0.0020791 

0.0006431 

0.0020849 

0.0008341 


A7--i»l-Node  Interface  With  Poisson  Ratio  of  0.2 


-0.4977104  0.2122934  0.0844444  0.0298734  0.0200340  0.0134027  0.0100929  0.0078050  0.0042497  0.0051307 

0.0042943  0.0034S44  0.0031490  0.0027424  0.0024101  0.0021344  0.0019030  0.0017044  0.0015380  0.0013924 

0.0012437  0.0011544  0.0010547  0.0009499  0.0008927  0.0008237  0.0007419  0.0007043  0.0004543  0.0004113 

0.0005707  0.0005344  0.0005021  0.0004737  0.0004494  0.0004300  0.0004144  0.0004123  0.0004213  0.0005044 

0.0007537 

-0.388252V  0.1741930  0.0719040  0.0213284  0.0142344  0.0092981  0.0047440  0.0051394  0.0040414  0.0032984 

0.0027377  0.0023117  0.0019794  0.0017130  0.0015003  0.0013234  0.0011740  0.0010514  0.0009431  0.0008539 

0.0 007747  0.0007054  0.0004449  0.0005913  0.0005439  0.0005014  0.0004440  0.0004303  0.0004000  0.0003731 

0.0003490  0.0003279  0.0003093  0.0002937  0.0002813  0.0002729  0.0002709  0.0002779  0.0003344  0.0003044 

-0.4024443  0.1711134  0.0484971  0.0190154  0.0125183  0.0079710  0.0057013  0.0042471  0.0033311  0.0024773 

0.0022024  0.0018437  0.0013701  0.0013324  0.0011773  0.0010340  0.0009131  0.0008154  0.0007309  0.0004586 

0.0005941  0.0005420  0.0004944  0.0004330  0.0004143  0.0003839  0.0003550  0.0003293  0.0003066  0.0002864 

0.0002687  0.0002534  0.0002407  0.0002307  0.0002241  0.0002230  0.0002294  0.0002779  0.0004213 

-0.4044387  0.1699114  0.0676011  0.0184101  0.0120501  0.0075963  0.0053937  0.0040096  0.0031121  0.0024887 

0.0020384  0.0017016  0.0014426  0.0012389  0.0010754  0.0009423  0.0008323  0.0007403  0.0006626  0.0005961 

0.0003391  0.0004897  0.0004467  0.0004090  0.0003759  0.0003466  0.0003209  0.0002781  0.0002780  0.0002606 

0.0002434  0.0002331  0.0002234  0.0002171  0.0002164  0.0002230  0.0002709  0.0004123 

-0.4032449  0.1493441  0.0472731  0.0180943  0.0117990  0.0073901  0.0052211  0.0038629  0.0029859  0.00237B7 

0.0019419  0.0016141  0.0013444  0.0011707  0.0010143  0.0008871  0.0007823  0.0006949  0.0006210  0.0005583 

0.0005046  0.0004580  0.0004177  0.0003824  0.0003516  0.0003246  0.0003009  0.0002801  0.0002620  0.0002467 

0.0002339  0.0002241  0.0002179  0.0002171  0.0002241  0.0002729  0.0004164 

-0.4034143  0.1490887  0.0470423  0.0179239  0.0114591  0.0072731  0.0051217  0.0037773  0.0029113  0.0023134 

0.0018841  0.0013447  0.0013204  0.0011294  0.0009770  0.0008534  0.0007517  0.0006670  0.0005957  0.0005333 

0.0004836  0.0004390  0.0004003  0.0003667  0.0003376  0.0003120  0. 0002899  0.0002706  0.0002543  0.0002409 

0.0002306  0.0002241  0.0002234  0.0002307  0.0002813  0.0004300 

-0.4038204  0.1689309  0.0669360  0.0178223  0.0115724  0.0071994  0.0050584  0.0037221  0.0028630  0.0022707 

0.0018461  0.0015309  0.0012901  0.0011021  0.0009323  0.0008310  0.0007314  0.0006486  0.0005791  0.0005201 

0.0004699  0.0004247  0.0003894  0.0003570  0.0003290  0.0003049  0.0002840  0.0002o44  0.0002519  0.0002409 

0.0002339  0.0002331  0.0002407  0.0002937  0.0004494 

-0.4039424  0.1688333  0.0448560  0.0177534  0.0115157  0.0071507  0.0050160  0.0036850  0.0028303  0.0022416 

0.0018201  0.0013074  0.0012493  0.0010833  0.0009333  0.0003137  0.0007176  0.0006341  0.0005479  0.0005101 

0.0004610  0.0004189  0.0003824  0.0003314  0.0003246  0.0003016  0.0002823  0.0002444  0.0002543  0.0002447 

0.0002456  0.0002534  0.0003093  0.0004737 

-0.6040204  0.1687694  0.0668026  0.0177101  0.0114767  0.0071169  0.0049863  0.0034589  0.0028071  0.0022210 

0.0018017  0.0014911  0.0012544  0.0010699  0.0009233  0.0008049  0.0007080  0.0006276  0.0005603  0.0005036 


A8— 41-Node  Interface  With  Poisson  Ratio  of  0.3  (Cont.) 


A7 


0 >0004533  0.0004141  0.0003780  0.0003487  0.0003231 
0.0002487  0.0003270  0.0005021 

-0.4040720  0.1487254  0.0447454  0.0174783  0.0114401 
0.0017884  0.0014701  0.0012437  0.0010404  0.0000140 
0.0004521  0.0004120  0.0003777  0.0003487  0.0003244 
0.0003400  0.0005344 

-0.4041004  0.1484047  0.0447380  0.0174553  0.0114200 
0.0017787  0.0014704  0.0012341  0.0010537  0.0000000 
0.0004511  0.0004120  0.0003780  0.0003514  0.0003200 
0.0003707 

-0.4041354  0.1484723  0.0447104  0.0174384  0.0114143 
0.0017717  0.0014444  0.0012307  0.0010401  0.0000051 
0.0004321  0.0004141  0.0003824  0.0003570  0.0003374 
-0.4041547  0.1484537  0.0447051  0.0174230  0.0114033 
0.0017447  0.0014401  0.0012271  0.0010441  0.0009020 
0.0C04S33  0.0004180  0.0003804  0.0003447  0.0003514 
0.4041400  0.1484434  0.0444044  0.0174144  0.0113050 
0.0017434  0.0014574  0.0012251  0.0010440  0.0000021 
0.0004410  0.0004247  0.0004003  0.0003824  0.0003750 
-0.4041707  0.1484341  0.0444844  0.0174004  0.0113880 
0.0017414  0.0014541  0.0012244  0.0010440  0.0000020 
0.0004400  0.0004300  0.0004177  0.0004000  0.0004143 
•0.4041874  0.1484273  0.0444804  0.0174043  0.0113844 
0.0017400  0.0014541  0.0012251  0.0010441  0.0000051 
0.0004834  0.0004380  0.0004447  0.0004330  0.0003430 
-0.4041034  0.1484221  0.0444741  0.0174007  0.0113814 
0.0017414  0.0014374  0.0012271  0.0010401  0.0000000 
0.0003044  0.0004807  0.0004044  0.0005013  0.0008027 
-0.4041070  0.1484187  0.0444733  0.0173084  0.0113800 
0.0017434  0.0014401  0.0012307  0.0010337  0.0000140 
0.0005301  0.0005420  0.0004440  0.0000400 
-0.4042007  0.1484144  0.0444714  0.0175073  0.0113704 
0.0017447  0.0014444  0.0012341  0.0010404  0. 0000233 
0.0003041  0.0007054  0.0010547 

-0.4042023  0.1484153  0.0444710  0.0175073  0.0113800 
0.0017717  0.0014704  0.0012437  0.0010400  0.0000333 
0.0007747  0.0011544 

-0.4042020  0.1484153  0.0444714  0.0175084  0.0113814 
0.0017787  0.0014701  0.0012544  0.0010833  0.0000523 
0.0012457 


0.0003014  0.0002840  0.0002704  0.0002420  0.0002404 

0.0070027  0.0u40651  0.0036401  0.0027004  0.0022041 
0.0007074  0.000701J  0.0006210  0.0005554  0.0004006 
0.0003040  0.0002800  0.0002801  0.0002780  0.0002864 

0.0070751  0.0040406  0.0036263  0.0027781  0.0021051 
0.0007023  0.0006070  0.0006183  0.0005527  0.0004076 
0.0003120  0.0003000  0.0002081  0.0003066  0.7003731 

0.0070420  0.0040380  0.0036161  0.0027601  0.0021873 
0.0007801  0.0006046  0.0006167  0.0005510  0.0004076 
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0.0007023  0.0007013  0.0006276  0.0005670  0.0005201 
0.0008237 
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0.0007074  0.0007080  0.0004341  0.0005701  0.0005353 
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0.0008040  0.0007174  0.0004484  0.0005057  0.0005583 
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0.0008310  0.0007517  0.0006040  0.0006626  0.0004584 

0.0070363  0.0040187  0.0036031  0.0027626  0.0021873 
0.0008534  0.0007823  0.0007403  0.0007300  0.0008530 


A9--*»1-Node  Interface  Node  With  Poisson  Ratio  of  0.3 


ooooooo  oooo 

Ui  UJ  UJ  UJ  UJ  Ui  UJ  UJ  LU  Ui  UJ 

no  o  o 

KN»>NTnnH  4<ONO 
Oo«'0'*N'«n  o.  >o  o 
*—  «  r*>  p*  ♦  «  «m  rv  in  ♦  p"  o 
«J**'M,CM*HCMI«IO  IP  O"  41  o 
«  o>  on  n  n«r  o  o  o 
IkNNnNMMO  OOOO 
n  N*<  «  Hrt  H  T+  *+  tH 

ooooooo  OOOO 


OOOO  O  1 
I  I  I  I  I 
111  UJ  UJ  UJ  Ui  I 

m  n  «i 

O  o  O  o  0»  I 
XNVMOw 
3  n  (M«>  n  mi 
W  If)  Oh  n  N 
•O  MN  ♦  0~ 

r**  o-  othi 
-O  JO  CM  CM  CM  I 


OOOO 
till 
Ui  Ul  UJ  UI 
•O  O  CM  00 

o*  o  o  ©>• 
o»  o  r*>  n 

HN  ♦  ». 
IT  N  H  fvj 
>0  H  (►  ffl 

oo  oo  rs  n 


ooooooo  oooo 

I  I  I  I  I  I  I  I  I  I  I 


oooo  o o 

UJ  UJ  UJ  UJ  UJ  i  ■  i 
•v  j  so-h  oo  o 
ecNoxon  no  u 

O  T  H00  CO  CM  O  o 

K  O*  O  H  P*»  ^  O  (\) 

U  in  NN  (V  (DO  «t_ 

«  <r  nooo  ho  v. 

U.  N  <ONO  OO  V 

(M  HH  H  HH  4J 

•  •»•••  C 

O  O  O  O  O  O  — 


OOOO  o  O 
I  I  I  I  II 
Ul  Ul  Ui  UJ  UJ  Ul 
in  «■*  'O  o  « 
O  O  O  O'  o  o 
z  n  «hm  cm o 
3«ooo  o>n 
coo>09r»-«  oo  cj 

O  &•  O'  H  H4 

♦  HO 

oo»wn  wn 
•  «  *  •  •  « 
oooo  o  o 

till  II 


o  o  o  o  o  o 

UJ  UJ  UJ  UJ  UJ  UJ 
N  N  ^  O 

KMOOOm  «o  4) 

O  H  O-ffl  H  O.  o  o 

>—  vo  *h  p~  o  ro  o  ro 

O  O-  CON  ♦  CO  O  V(_ 

«  rt  MO  o  H  O  u 

U.  M  4N  o  O  O  4) 

MHHH  HH  4J 

o  o  o  o  o  o  — 


©oooo© 
I  I  I  I  II 
UJ  Ul  Ul  UJ  Ul  UJ 

T  H  pi  M  »•  M 

0  0  0  0  0-0. 

s  n  o  o  o  o>  r  j 

□  M  M  <■  N  CM  1*1 

in  O  ON  M  NO 

»  oh  Ml  N  N 

N  O  OJ  03 

N  n  NN 

•  ••••♦ 
o  o  o  ©  ©  © 
I  I  I  I  I  I 


v*v  v  v\* 

'j*  *  j  *  ja  "  (  •  *  « 


SUM  FACTOR  SUM  FACTOR 

-0.63916201E-01  0.4589275AE  01  -0.70285104E-01  0.45953288E  01 


W  111 


rH 

«-l  H 

rH 

H  H 

4-4  »4 

H 

rH 

H 

rH 

—4  4-4 

^H 

rH 

r-4 

4-4 

4-H 

^H 

rH 

r-4 

4-4 

© 

o  o 

© 

©  © 

o 

©  © 

o 

O 

o 

o 

o  o 

© 

o 

O 

© 

© 

© 

© 

© 

O 

© 

© 

© 

© 

© 

o 

© 

Ul 

Ul  Ul 

Ui 

Ul  Ul 

Ui 

Ul  Ul 

UJ 

UJ 

UJ 

UJ 

Ul  UJ 

UJ 

UJ 

UJ 

Ul 

UJ 

UJ 

UJ 

UJ 

UJ 

Ui 

Ul 

UJ 

UJ 

Ul 

Ul 

Ul 

rv 

o-  O' 

-o 

nr  cm 

00 

O-  IV 

o 

o- 

in 

hp 

CM  00 

UO 

Hp 

PO 

rH 

O' 

P0 

UO 

CM 

'O 

© 

O' 

O' 

- o 

ro 

CM 

© 

0£ 

m3 

ri  -o 

rv  wo 

CM 

4-4  -O 

-o 

«r 

03 

-o 

cm  po 

© 

■V 

4*H 

CM 

00 

CJ 

wo 

Pv 

P' 

rv 

00 

CM 

-o 

© 

CM 

o 

O 

ro 

-h  fO 

o 

o*  -o 

PV 

Pv  o 

ro 

o» 

O' 

wo 

po  -o 

o 

GO 

4-4 

O' 

rH 

03 

’O 

O' 

rH 

PO 

CM 

"O 

rv 

O' 

© 

y 

H* 

00 

o-  rH 

-o 

M3  >0 

>o 

oo  4- 

00 

-o 

"O 

o 

■o  4 

O' 

P0 

UO 

GO 

ro 

UO 

HP 

fM 

in 

UO 

IV 

© 

UO 

O' 

O’ 

© 

(0 

CJ 

o 

«-•  'O 

-© 

ro  -o 

CM 

sO  o 

rv 

wo 

rH 

O’ 

4-H  CO 

CO 

rv 

© 

rH 

IV 

PO 

PO 

UO 

Ctk  CM 

» 

rH 

o 

© 

it- 

<c 

in 

M  o 

CM 

00  rH 

O' 

a-  cm 

m 

o 

nO 

CM 

O'  'O 

ro 

rH 

O' 

00 

'O 

in 

o- 

PO 

CM 

rH 

rH 

© 

o 

© 

© 

© 

u. 

HP 

in  in 

ri 

00  Ps 

UO 

*r  «r 

ro 

ro 

CM 

CM 

4-4  4-4 

rH 

rH 

© 

© 

© 

© 

© 

© 

© 

© 

© 

o 

o 

o 

o 

© 

H) 

'O 

PO  (M 

CM 

H 

4-4 

r4  4-4 

»4 

4-4 

4-H 

•H  1"H 

^4 

4-4 

^H 

4-4 

4-4 

4^ 

rH 

rH 

rH 

rH 

4^4 

rH 

^H 

4-H 

4-H 

rH 

4-f 

© 

©  © 

O 

©  © 

o 

o  o 

o 

© 

o 

o 

o  © 

o 

o 

o 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

o 

© 

© 

c 

iH 

*-«  rH 

*■4 

»H 

4-4  4-4 

rH 

rH 

rH 

rH 

4-1  -I 

t4 

4* 

4-4 

rH 

r-4 

rH 

rH 

rH 

0) 

© 

©  © 

o 

©  © 

o 

o  o 

© 

o 

© 

o 

©  © 

o 

©  © 

© 

o  o 

o 

© 

o 

o 

© 

© 

© 

© 

© 

© 

o 

Ul 

Ul  UJ 

UJ 

UJ  UJ 

UJ 

ui  ui 

Ul 

UJ 

Ui 

Ul 

UJ  Ui 

Ul 

Ui  Ul 

Ul 

Ul  Ul 

Ul 

Ul 

ill 

ill 

Ul 

UJ 

ui 

Ul 

ui 

1 

UJ 

z 

1 

rH 

m 

N 

rH  Ul 

o 

in  o» 

cm 

CM 

pv 

ro 

r-4  CM 

00 

'C  CJ 

Pv 

©  © 

rv 

UO 

UO 

Hp 

00 

© 

O' 

«T 

wo 

CM 

© 

o  o 

O'* 

o  o 

o 

©  O' 

o 

O 

O'* 

o 

O  o 

O' 

O'  © 

O' 

O'  © 

O' 

O' 

O' 

O' 

O' 

© 

O' 

O' 

O' 

© 

x: 

»-4 

CM  CM 

in 

©  IV 

o 

fM  MT 

*r 

o 

M3 

in 

CM  O* 

O' 

MD  O' 

CM 

©  UO 

UO 

o 

rH 

>a 

00 

'T 

CO 

rv 

o 

©  in 

IV.  oo 

— «  ro 

CM 

O'  in 

rv 

o- 

O' 

UO 

4-4  PO 

hp 

O 

PT 

H 

M3  UO 

>o 

fV 

M3 

M 

nr 

CM 

0“ 

© 

o 

tn 

o 

<o  -o 

00 

m  o. 

P3 

4  CM 

00 

00 

O' 

O' 

-o  co 

UO 

>o 

rH 

O' 

O'  CM 

p* 

o* 

PO 

<P 

-o 

© 

UO 

CJ 

© 

O* 

CM 

*<r 

*-<  in 

o* 

CM  O' 

ro 

©  o 

4-4 

CO 

PO 

O'  UO 

CM 

O' 

?V 

Hp 

CM 

rH 

O' 

03 

rv 

'O 

wo 

WO 

*r 

•r 

ro 

■o 

fv  i-o 

00 

ill  CM 

n 

O  O' 

03 

r-v 

-O 

M3 

in  ui 

toO 

■ht 

«f 

4T 

O' 

Hp 

ro 

ro 

PO 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

GQ 

00 

<r  ro 

CM 

cm  rM 

CM 

CM  *-< 

4-4 

f-4 

4— < 

+  4  4-H 

rH 

—H 

r-4 

rH 

4-H 

rH 

rH 

r4 

rH 

rH 

rH 

rH 

^H 

4-H 

— < 

©  ©  ©  o  o  o 

lit  III 


©  o  o 

I  I  I 


OOOOOOOOOOOOOOOOOOOOO 

I  I  I  I  '  I  I  I  !  I  !  I  I  I  i  I  i  I  I  ; 


©000000000000000000000000000000 


QC. 

UJ 

f% 

UJ  UJ 
CM  U0 

ro  c* 

UJ 

PO 

CM 

Ul  Ul 
CM  CM 
00  IV 

Ul 

O' 

CO 

UJ  Ui 

o  o 

rH  ro 

UJ 

© 

Ui  UJ 
r  t'j 
CM  00 

Ul 

© 

© 

Ul  ui 
ro  © 
PO  CD 

UJ 

in 

O' 

UJ  UJ 

r-4  Pv 
©  © 

UJ 

P' 

Ul  Ul 

ro  © 
in  © 

ui 

PO 

in 

UJ 

© 

CM 

UJ 

CM 

in 

Ul 

CM 

M3 

Ul  UJ 

in  © 

©  4-H 

UJ 

© 

UJ 

PO 

«r 

UJ 

CM 

in 

UJ 

© 

© 

O 

rH 

CO  CM 

O' 

CM  rv 

CM  *“« 

CM 

rv  cm 

T  T 

© 

rv  ir 

rv 

CM  O 

© 

© 

© 

CM 

00  CM 

ro 

CM 

© 

H 

CP*  rH 

CO 

CM  © 

CO 

a  in 

PO 

P' 

UO 

O'  '■o 

O' 

ro  -«r 

© 

cm  in 

■«r 

Pv 

in 

in 

fV  rH 

in 

© 

in 

o 

0) 

C3 

ro 

O'  rH 

ro  cm 

'O 

O'  CM 

O' 

O  cm 

<T 

rH  ro 

© 

(V  00 

rH 

rv  <r 

ro 

ro 

UO 

00 

CM  © 

CM 

o 

© 

o 

<s  IV 

a  cm 

r 

O  CM 

O' 

O'  CM 

UO 

o  >o 

CM 

O'  "O 

ro 

rH  O' 

© 

'O  wo 

ro 

CM 

rH 

rH  O 

© 

© 

© 

© 

P0 

u. 

CM 

in  m 

rH 

00  Pv 

UO 

r  r 

PO 

PO  CM 

CM 

rH 

rH 

rH  © 

© 

©  o 

© 

© 

© 

o 

O  o 

© 

© 

© 

© 

4- 

>o 

ro  cm 

CM 

4-H  4-4 

rH 

rH  4-4 

4-H 

*4  rH 

rH 

rH  tH 

4-H 

rH  r-4 

rH 

4-H  rH 

rH 

rH 

rH 

^H 

4-H  rH 

4-4 

4-H 

rH 

4-H 

U 

o 

©  © 

o 

©  © 

© 

©  © 

o 

©  o 

o 

©  © 

© 

©  © 

© 

o  © 

© 

© 

© 

© 

©  © 

© 

© 

© 

© 

Q) 

4-1 

c 

rH 

r-4  rH 

—4 

rH  rH 

4-H 

rH  rH 

rH 

4-4 

4^ 

r-»  H 

rH 

4-H  4-H 

rH 

rH  rH 

4-H 

rH 

rH 

rH  rH 

rH 

4-H 

4*H 

4-H 

o 

©  o 

© 

©o 

© 

©  © 

© 

?? 

O 

?? 

© 

©  © 

© 

©  o 

© 

1 

Ui 

o 

o 

© 

o  © 

© 

© 

© 

© 

-g 

Ul 

UJ  UJ 

UJ 

U!  Ul 

UJ 

Ul  Ul 

UJ 

ui  UJ 

Ul 

Ul  Ul 

UJ 

Ul  Ul 

UJ 

Ul  UJ 

Ul 

Ul 

UJ 

UJ  UJ 

Ul 

Ul 

til 

UJ 

o 

in  o 

CD 

4-h  rr 

o 

ro  O' 

o 

O'  O' 

rH 

m  iv 

PO 

UO  © 

•33 

©  wo 

rH 

NO 

rv 

CM 

rH  >Q 

rH 

© 

CM 

03 

z 

o 

©  © 

O' 

©  © 

© 

O  O' 

O' 

O'  o 

© 

O'  O' 

O' 

O'  © 

O' 

O'  O' 

© 

O' 

O' 

© 

©  O' 

© 

© 

© 

O' 

1 

s 

-4T 

©  P' 

in 

ro 

O' 

©  *“• 

© 

n  m 

CM 

©  rv 

O' 

rH  Pv 

« 

uo  ro 

O' 

ro 

<r 

rv 

PO  V© 

rv 

rv 

CM 

*T 

3 

CM 

O'  © 

ir 

4-H  rH 

CM 

oo  O' 

CM 

CM  O 

O' 

iio  rv 

O' 

in  O' 

O' 

©  CM 

4-H 

Pv 

IV 

4-H 

O'  oo 

O' 

CJ 

M3 

O 

\0 

1 

co  n 

CO  CO 

O' 

'O  4-H 

CO 

'O 

O-  4-4 

O' 

r  ro 

o 

ro  cm 

«c 

O'  uo 

ro 

CM 

ro 

-o 

O' 

© 

00 

M3 

>o 

rt 

O'  P0 

rj 

U0  Pv 

ro 

cm 

p' 

rv 

CM 

O'  4 

ro 

4-4  O' 

pv 

m  -r 

ro 

CM 

4-H 

© 

O'  O' 

O' 

03 

C3 

CD 

00 

CD  pv 

PO 

©  © 

p' 

M3  UO 

•<  ro 

ro 

CM  CM 

CM 

CM 

rH 

4-H  rH 

4-H 

rH 

rH 

rH 

o  o 

© 

o 

© 

© 

*— 

'O 

ro  cm 

Cw| 

CM  rH 

rH 

•M  4-4 

rH 

4-P  rH 

4-4 

4-H  rH 

•H 

TH  4-4 

rv 

4-4  rH 

rH 

— * 

4-4 

rH 

4-4  r«0 

tH 

4-4 

4-H 

4-H 

CD 

o 

©  o 

© 

©  © 

© 

©  © 

© 

©  o 

o 

o  © 

© 

©  O 

© 

©  © 

© 

© 

© 

© 

©  o 

o 

© 

O 

© 

I  I  I 


© 

©  o 

© 

©  O 

© 

©  © 

o 

© 

© 

o 

©  © 

O 

© 

© 

o 

o 

o 

© 

© 

o 

o 

o 

© 

© 

©  © 

o 

UJ 

UJ  UJ 

Ul 

UJ  ui 

UJ 

Ul  Ui 

UJ 

Ul 

UJ 

UJ 

ul  ul 

UJ 

Ul 

UJ 

UJ 

UJ 

Ul 

UJ 

UJ 

UJ 

Ul 

UJ 

UJ 

Ui 

UJ  UJ 

UJ 

o 

03  O' 

-o 

O'  in 

p^ 

cm 

rH 

o- 

O' 

CM 

rH  ro 

CM 

CM 

© 

PO 

ro 

rv 

rH 

•** 

CM 

VI 

M3 

'O 

-o 

O'  PO 

© 

oc 

00 

in  co 

rH 

o-  ro 

'O 

cm  -r 

CM 

© 

rH 

o  a 

00 

M3 

O' 

CM 

40 

in 

Pv 

CM 

rv 

CM 

rH 

rv  o 

© 

o 

Ul 

PO  >o 

© 

rH  CM 

ui 

CM  CM 

rH 

ro 

ro 

CM 

«r  in 

'O 

© 

in 

Ul 

rH 

© 

4-4 

IT 

in 

© 

in 

PV 

PO  PO 

o 

41 

ro 

MJ  CM 

ro 

^  P' 

ro 

rH  -O 

P0 

rv 

ro 

rv  r 

rv 

CM 

CM 

'O 

rH 

ro 

CM 

o 

ro 

o* 

'O 

© 

O'  HP 

o 

o 

CJ 

CM  HP 

O' 

-r  »*> 

v 

©CM 

O' 

CM 

HT 

r  P0 

CO 

PV 

03 

P' 

ro 

ro 

too 

00 

CM 

CO 

o* 

4*  O 

© 

10 

« 

-o 

GO  CM 

o- 

o*  CM 

O' 

©  CM 

U! 

© 

'O 

CM 

O'  ^ 

PO 

^H 

O' 

03 

M3 

in 

4t 

ro 

CM 

rH 

4*H 

© 

© 

O  © 

© 

u- 

u. 

CM 

in  Ul 

rH 

00  iv 

ui 

in  <r 

ro 

ro 

CM 

CJ 

rH  ^H 

^H 

4-4 

© 

© 

© 

© 

© 

© 

© 

© 

o 

© 

© 

o  © 

o 

u 

PO  CM 

CM 

rH  »H 

rH 

4-4  rH 

rH 

rH 

^H 

rH 

rH  4-H 

rH 

rH 

4-H 

4-H 

4^ 

4-H 

^H 

rH 

rH 

rH 

4-H 

4-H 

rH 

rH  ^4 

rH 

41 

© 

©  © 

© 

©  © 

© 

©  © 

o 

©  © 

© 

o  o 

© 

© 

© 

© 

© 

© 

© 

o 

© 

© 

© 

© 

© 

©  © 

© 

c 

rH  4-H 

4—4 

rH  rH 

—4 

H  4-H 

rH 

rH 

rH 

rH 

rH  4-4 

rH 

rH 

4-H 

4-H 

rH 

4-H 

4-H 

CM 

CM 

CM 

CM 

CM 

CM  CM 

CM 

41 

? 

©  © 

1  1 

© 

1 

©  © 

1  1 

© 

1 

©  © 

1  1 

© 

1 

©  © 

1  t 

© 

1 

©  © 

f  1 

© 

T 

© 

1 

© 

1 

© 

1 

O 

1 

© 

t 

© 

1 

© 

1 

© 

l 

© 

1 

© 

1 

© 

1 

© 

1 

©  © 

1  1 

© 

i 

X) 

o 

Ul 

Ul  Ui 

UJ 

ui  uj 

Ul 

Ui  Ui 

UJ 

Ui  ui 

UJ 

UJ  ui 

Ul 

Ui 

Ul 

Ui 

Ui 

UJ 

Ul 

Ul 

Ui 

UJ 

Ul 

Ul 

UJ 

Ul  Ul 

Ui 

z 

'O 

rH  © 

ro 

P0  O' 

O' 

H  CM 

o 

O' 

0“ 

O' 

po  in 

rH 

CM 

4-i 

in 

© 

4-H 

© 

© 

Ul 

Pv 

ro 

O' 

©  O' 

© 

1 

© 

o  © 

O' 

O  O' 

o* 

o© 

O' 

o 

© 

O' 

O  CP- 

© 

©  O' 

© 

O' 

© 

© 

© 

CM 

ro 

CO 

rH 

CO 

rv  O' 

in 

2C 

Ul 

Pv  rH 

o 

O'  © 

«> 

4-»  H) 

ro 

p*- 

© 

P3 

05  © 

O* 

Cm 

Hf 

-o 

CO 

© 

i> 

>o 

© 

© 

O' 

o 

O' 

O'  O' 

© 

vO 

"Zi 

O' 

i n  a> 

CM 

rv  po 

rH 

^c  CO 

r4 

03 

n 

M3 

rv  CM 

n 

^4 

< 

<r 

CM 

00 

Ov 

'O 

GO 

o 

O' 

PO 

© 

nr  rv 

'O 

1 

CD 

ro 

CM  © 

in 

vO  © 

Pv 

CM  -C 

Ul 

PO 

© 

ro 

rH  «T 

© 

O' 

© 

HT 

© 

rv 

vO 

Pv 

Ul 

CD 

rH 

rv 

CM 

Pv  CM 

CO 

1 

© 

O'  "O 

o* 

«r  oo 

MO 

CO 

Cl 

r> 

ro 

O' 

«o  PO 

«— 

03 

Pv 

m 

«r 

CM 

4—4 

© 

O' 

CM 

rv 

CM 

O' 

>e  wi 

Hp 

o 

4-H 

© 

00  'O 

in 

o*  m 

Pi 

PM 

CM 

rH 

rH  • 

rH 

© 

© 

© 

© 

© 

© 

© 

CT' 

O' 

CO 

00 

pv 

<v  rv 

fV 

vO 

ro  cm 

CM 

r  rH 

— < 

H  n 

rH 

rH 

4-H 

4—  4-4 

s-H 

r  4 

rH 

— 

4-4 

rH 

4-H 

•> 

O' 

O' 

O' 

O' 

O'  O' 

O' 

CD 

o 

1 

©  © 

»  i 

© 

* 

©  © 
t  1 

© 

i 

©  o 

!  t 

»-V. 

mT 

© 

\ 

O 

1 

©  © 

;  1 

© 

© 

« 

© 

: 

© 

© 

© 

1 

© 

l 

© 

o 

1 

©  o  © 

1  1  ! 

o 

1 

o  © 

!  1 

© 

1 

B5 


'  A  *  ^  *  «  »  . 


iJUULfcfc 


-•  *g-f 


APPENDIX  C 


CENTERLINE  ROWS  FOR  61-NODE  INTERFACE 


APPENDIX  F 

41-NODE  BEAM  ON  FOUNDATION  WITH  POISSON'S  RATIO 
OF  0.3  AND  MODULUS  OF  ELASTICITY 
EQUAL  TO  7000  PS  I 


APPENDIX  E 


31 -NODE  BEAM  ON  FOUNDATION  WITH  POISSON'S  RATIO 
OF  0.3  AND  MODULUS  OF  ELASTICITY 


EQUAL  TO  7000  PS  I 


APPENDIX  D 


21-NODE  BEAM  ON  FOUNDATION  WITH  POISSON'S  RATIO 
OF  0.3  AND  MODULUS  OF  ELASTICITY 


EQUAL  TO  7000  PS  I 


♦56404965E  00  0.16103446E  00  0.58458386E-01  0.15996674E-01  0. 10742534E-01  0.64344953E-02  0.45185527E-02 
.  32922967E-02  0.25169003E-02  0, 19836183E-02  0.16040107E-02  0. 13236628E-02  0.11109951E-02  0.94591B98E-03 
.81534584E-03  0.71041205E-03  0.&2495845E-03  0.55460253E-03  0.49616803E-03  0.44732110E-03  0.40633504E-03 
.37193824E-03  0.34321996E-03  0.31958583E-03  0. 30075588E-03  0.28688510E-03  0.27843818E-03  0.27888179E-03 
.28978657E-03  0.35464108E-03  0 .55824179E-03 
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APPENDIX  G 

21- NODE  BEAM  ON  FOUNDATION  WITH  POISSON'S  RATIO 
OF  0.2  AND  MODULUS  OF  ELASTICITY 
EQUAL  TO  7000  PS  I 


APPENDIX  H 

31-NODE  BEAM  OF  FOUNDATION  WITH  POISSON'S  RATIO 
OF  0.2  AND  MODULUS  OF  ELASTICITY 


EQUAL  TO  7000  PS  I 
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APPENDIX  I 


41-NODE  BEAM  ON  FOUNDATION  WITH  POISSON'S  RATIO 


OF  0.2  AND  MODULUS  OF  ELASTICITY 


EQUAL  TO  7000  PS1 
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APPENDIX  J 


SEQUENTIAL  SEARCH  METHOD  BY  NELDER  AND  MEAD 


The  sequential  search  method  was  used  to  find  the  coefficients  for 
given  functions  to  fit  a  row  of  the  foundation  stiffness  matrix  in  a 
least  squares  sense.  This  sequential  search  method  is  an  optimization 
process  known  as  the  flexible  simplex  search  developed  by  Nelder  and  Mead 
(23,  2k). 

This  procedure  considers  an  objective  function,  f,  defined  through¬ 
out  an  N-dimensional  space  with  coordinates  (X, ,  X.,  X,,  .  .  .,X  ).  The 

I  l  i  n 

search  begins  with  an  initial  point  and  generates  N  more  points  with  a 
total  of  N  +  1  corner  points  (called  a  simplex).  Each  corner  point  defin 
ed  in  the  N-dimensional  space  is  systematically  replaced  by  new  ones  in 
order  to  make  the  objective  function  produce  a  minimum  value.  As  the 
method  produces  points  which  are  producing  smaller  objective  function 
evaluations,  the  simplex  becomes  smaller  until  the  desired  minimum  is 
obta i ned . 

The  interactive  procedure  of  replacing  the  corner  points  uses  three 
moves  called  reflection,  expansion,  and  contraction.  For  each  trial 
point  which  is  a  corner  to  the  simplex,  the  objective  function  is  evalu¬ 
ated.  The  corners  are  then  numbered  Z  .  .  .  Z  ,  where  the  smallest  func 

o  n 

tion  value  corresponds  to  Fq  and  the  largest  function  value  corresponds 

to  F  .  The  centroid,  Z,  is  calculated  from  pointsZ,...,Z  A  re- 

n  o  n- 1 

fleeted  point,  Zr>  is  calculated  as  follows: 

Z  =*  Z  +  A(Z  -  Z  ) 
r  n 

where  A  is  a  reflection  coefficient  (positive  number).  The  objective 
function  evaluation  of  this  point  is  labeled  Fr: 


The  magnitude  of  F^  is  examined  to  determine  further  steps  accord¬ 
ing  to  the  following  cases: 

Case  1 :  F  <  F  <  F  . 

o  r  n-1 

Case  2:  F  <  F 
r  o 

Case  3:  F  .  <  F  . 

n  - 1  r 

For  Case  I,  replace  Z with  point  Zr  and  begin  iteration.  For  Case 
2,  a  better  point  has  been  found  and  an  attempt  is  made  to  expand  fur¬ 
ther  along  the  same  direction.  The  expanded  point,  Z  ,  is  calculated  as 

6 

fol lows: 

Z  =  Z  +  B{Z  -  Z) 
e  n 

where  B  is  an  expansion  coefficient  (positive  number).  The  objective  is 

evaluated  for  Z  to  determine  F  .  If  F  <  F  ,  the  expansion  is  accepted, 

e  e  e  o 

Then  Z  is  replaced  with  Z  and  iteration  is  begun.  If  F  >  F  ,  the  ex- 
n  e  e  o 

pansion  is  rejected  and  Z n  is  replaced  with  Z^.  Then  a  new  iteration  is 
started. 

For  Case  3,  the  simplex  is  too  large  and  a  contraction  point,  Z c» 
is  calculated  as  follows: 

For  F  <  F  :  Z  =  Z  +  C(Z  -  Z) 
n  r  c  n 

For  F  >  F  :  Z  =  Z  +  C(Z  -  Z) 
n  r  c  r 

where  C  are  contraction  coefficients.  The  objective  function  is  then 

evaluated  to  determine  F.  IfF  <F  <F  or  F  <F  <F,  the  contrac 

c  c  n  r  c  r  n 

tion  is  completed  by  replacing  Zn  with  Z and  starting  iteration.  Other 
wise,  contraction  is  made  by  moving  each  corner  point  halfway  toward  the 


point  with  the  lowest  function  evaluation.  Then  a  new  iteration  begins. 


APPENDIX  H 

41 -NODE  BEAM  ON  FOUNDATION  GENERATED  FROM 
COEFFICIENTS  WITH  POISSON'S  RATIO  OF 
0.3  AND  MODULUS  OF  ELASTICITY 
EQUAL  TO  7000  PS  I 
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APPENDIX  L 

31-NODE  BEAM  ON  FOUNDATION  GENERATED  FROM 
COEFFICIENTS  WITH  POISSON'S  RATIO  OF 
0.3  AND  MODULUS  OF  ELASTICITY 
EQUAL  TO  7000  PS  I 
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APPENDIX  K 

21-NODE  BEAM  ON  FOUNDATION  GENERATED  FROM 
COEFFICIENTS  WITH  POISSON'S  RATIO  OF 
0.3  AND  MODULUS  OF  ELASTICITY 
EQUAL  TO  7000  PS  I 


A 


w 

x 

y 


coefficients  for  exponential  equation 
coefficients  for  exponential  equation 
coefficients  for  exponential  equation 
coefficients  for  exponential  equation 
second  soil  parameter 
approximate  interface  stiffness  matrix 
force  applied  to  beam-column 
modulus  of  subgrade  reaction 
beam-column  stiffness  matrix 

stiffness  term  for  nodes  on  interface  between  FEM  foundation  and 
beam-column  beam 

interface  stiffness  matrix  which  contains  all  characteristics  of 
the  two-dimensional  foundation 

stiffness  terms  coupling  the  interface  to  the  remainder  of  the 
foundation 

stiffness  terms  coupling  the  interface  to  the  remainder  of  the 
foundation 

stiffness  terms  of  remainder  of  foundation 
surface  pressure 

beam-column  (interface  nodes)  displacements 
.displacement  of  remainder  of  foundation 

surface  displacements 

number  of  nodes  from  major  diagonal 
magnitude  of  interface  stiffness  term 
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